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ABSTRACT 


^ Time  series  measurements  of  shallow-water  waves  and  fluid-mud  density  variations, 
taken  simultaneously  with  tide  elevation  and  suspended  sediment  data,  indicate  that 
wave/fluid-mud  interactions  in  the  nearshore  may  be  largely  responsible  for  the  present-day 
accumulation  of  fine-grained  sediment  on  the  open,  unprotected  coast  of  northeastern  South 
America.  Results  of  field  experiments  conducted  along  the  central  Surinam  coast  show  that 
accumulations  of  fluid  mud,  which  can  be  found  up  to  1.5  m thick  on  western  flanks  of 
migrating  shore-attached  mudbanks,  affect  incoming  swell  by  changing  their  form  from 
sinusoidal  to  solitary-like  and  by  preventing  wave  breaking  except  for  occasional  spilling.  As 
long-period  swell  (T  = 12-16  sec)  propagates  over  shallow  banks  of  fluid  mud,  the  wave  height 
to  water  depth  ratio  remains  nearly  constant  at  0.23;  the  steady  decrease  in  wave  height  with 
shoaling  water  depth  indicates  that  substantial  amounts  of  wave  energy  are  lost  to  a fluid-mud 
bottom  even  though  breaking  does  not  occur.  — — 

Evaluation  of  several  mechanisms  for  the  dissipation  of  wave  energy  reveals  that 
attenuation  of  high-frequency  waves  (T  < 3 sec)  may  be  due  to  internal  friction  resulting  from 
the  extraordinarily  high  kinematic  viscosity  of  muddy  coastal  waters.  Additional  energy  is  lost 
by  viscous  shear  within  a bottom  boundary  layer.  Agreement  with  Keulegan's  (1948)  theory  for 
energy  dissipation  for  solitary  waves  on  a smooth  bottom  is  demonstrated  when  near-bottom 
boundary  layer  viscosities  of  0.1  to  1.0  stoke  are  taken  on  a fluid-mud  bottom.  Additionally, 
unknown  but  large  amounts  of  energy  are  lost  as  waves  suspend  soft,  fluid  mud. 

A specially  designed  wave/fluid-mud  pressure-sensing  system,  used  successfully  to 
obtain  mean  and  fluctuating  components  of  near-bottom  or  subbottom  mud  density  variations 
during  the  passage  of  surface  water  waves,  indicates  that  fluid-mud  density  fluctuates  by  0.025 
to  0.05  gm/cm^  on  a wave-by-wave  basis.  As  each  wave  crest  passes  the  sensing  instrument,  a 
cloud  of  sediment  with  particle  concentration  of  approximately  50,000  mg/1  is  suspended 
instantly;  rapid  settling  that  follows  is  explained  by  the  formation  of  a collapsing  interface. 
Wave-by-wave  sediment  suspension  is  superimposed  on  a longer  term  suspension  and  deposition 
process  that  is  related  to  stage  of  the  tide.  In  water  1-5  m deep,  suspended  sediment  concen- 
trations increase  beginning  near  time  of  mid-tide,  reach  a maximum  at  low  water,  then 
decrease  steadily  up  to  the  following  high  water.  Surface  suspensate  concentrations  may 
exceed  3500  .ng/1  during  this  exchange  process  as  0. 1-0.8  m of  fluid  mud  is  suspended  in  a 3- 
hour  period  of  time  near  low  water. 


Utilizing^he  concept  of  wave-associated  currents,  it  has  been  shown,  in  taking  angle  of 
wave  approach,  average  observed  concentration  of  sediment  in  suspension,  and  net  drift  as 
given  by  solitary  wave  theory,  that  volume  transports  oT 3 to  70  x 10^  mVyrTan  be  explained 
by  waves  aloneOvJhese  rates  are  perhaps  10-100  times  higher  than  volume  transport  rates 
along  typical  sandyVnvironments;  high  volume  transport  rates  found  by  taking  this  approach  do 
not  require  breaking\waves  and  the  concept  of  radiation  stress  and  a nearshore  circulation  cell, 
as  normally  applied  4o  sandy  coasts.^  If  a hypothesis  of  mud  transport  by  solitary  waves  is 
accepted,  then  the  continuous  shorewclrd  transport  of  suspended  fluid  mud,  combined  with  the 
high  rate  of  wave  energy  dissipation,  explains  in  part  the  ability  of  muddy  coasts  to  protect 
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INTRODUCTION 


A view  traditionally  held  by  geologists  is  that  mud  accumulates  only  in  quiet,  shallow 
basins;  yet,  the  longest  uninterrupted  mud  coast  in  the  world  spans  1600  km  between  the 
Amazon  and  Orinoco  Rivers  along  the  open,  unprotected  coast  of  South  America,  where 
deepwater  wave  statistics  indicate  high  wave-energy  levels.  Many  of  the  physical  and 
geological  processes  that  are  important  in  maintaining  this  muddy  environment  take  place  in 
the  nearshore  region,  where  large  accumulations  of  thixotropic  gel,  referred  to  as  fluid  mud  or 
slingmud,  are  present.  Recent  studies  have  shown  that  the  transportation  of  much  of  the 
sediment  along  this  coast  results  from  the  northwestward  migration  of  large  shoref ace- 
attached  mudbanks,  parts  of  which  are  composed  of  fluid  mud  (Delft  Hydraulics  Laboratory; 
1962;  Allersma,  1968;  NEDECO,  1968).  Further,  the  distribution  of  nearshore  wave  energy 
along  the  northeast  coast  of  South  America  is  controlled  largely  by  these  mudbanks,  since 
accumulations  of  fluid  mud  have  an  attenuating  effect  on  waves.  Many  of  the  unanswered 
questions  about  this  muddy-coast  environment  concern  the  details  of  suspension,  transport,  and 
accumulation  of  fluid  mud  and  the  interaction  between  fluid  mud  and  incoming  waves.  So  that 
this  interesting  environment  might  be  better  understood,  this  investigation  was  undertaken  in 
order  to  examine  the  process-form  interactions  between  waves,  fluid  mud,  and  suspended 
sediment  in  the  nearshore  region. 

Comparison  among  hydrologic  and  geologic  surveys  dated  as  early  as  1670  indicated  that 
large  banks  of  fluid  mud,  originating  near  the  mouth  of  the  Amazon  River,  migrated  to  the 
northwest  at  an  average  rate  of  1.5  km/yr  (Allersma,  1968).  Extensive  field  investigations  in 
British  Guiana  (Delft  Hydraulics  Laboratory,  1962)  and  Surinam  (NEDECO,  1968)  affirmed  the 
fact  that  these  mudbanks,  although  exposed  to  wave  forces  from  the  Atlantic  Ocean,  were  able 
to  maintain  their  identity  for  decades.  These  studies  further  established,  from  visual 
observations,  that  waves  were  attenuated  rapidly  when  moving  over  mudbanks,  often  to  the 
point  of  disapp>earing  before  reaching  shore.  Although  geomorphology  of  mudbanks, 
characteristics  of  coastal  sediments,  and  general  coastal  oceanography,  such  as  open  shelf 
circulation,  have  been  discussed  in  the  literature  (Zonneveld,  195'*;  Vann,  1959;  Delft 
Hydraulics  Laboratory,  1962;  Eisma,  1967;  Allersma,  1968;  NEDECO,  1968;  Gibbs,  1975,  1976), 
to  date  there  is  little  published  data  on  wave  statistics  and  wave/fluid-mud  interactions  in  the 
nearshore  region  of  northeastern  South  America. 

As  early  as  1968,  the  need  for  such  wave  data,  other  than  from  visual  estimation,  was 
noted  by  Allersma: 


...it  is  felt  that  one  or  more  series  of  observations  with  instru- 
ments, combined  with  a thorough  analysis  of  the  records,  might  reveal 
some  facts  that  cannot  be  obtained  from  visual  observations. 

Specifically  what  is  needed  are  continuous  time  series  measurements  of  wave  height,  period, 
form,  and  angle  of  approach.  Each  time  series  should  be  of  sufficient  length  to  obtain  the 
nearshore  wave  characteristics.  Then,  to  answer  specific  questions  concerning  mudbank 
formation,  wave  attenuation  and  details  of  fluid-mud  transport,  additional  data  on  nearshore 
currents,  tide  elevation,  bottom  changes,  and  suspended-sediment  variations  must  be  collected 
and  analyzed. 
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In  view  ol  the  above,  the  specific  objectives  of  this  study  were:  1)  to  describe  waves 
moving  over  a mobile,  fluid-mud  bottom  by  using  electronic  sensors  and  to  examine  the  effect 
of  this  type  of  bottom  and  fluid-mud  suspension  on  wave  form  and  2)  to  assess  the  role  of 
waves  in  suspending  and  transporting  fluid  mud.  These  objectives  were  accomplished  through 
field  experiments  conducted  in  four  different  hydrologic  settings  along  the  central  Surinam 
coast.  Much  of  the  data  on  wave/fluid-mud  interactions  was  collected  by  utilizing  a 
wave/fluid-mud  gage  that  was  designed  specifically  for  muddy  environments.  Time  series 
records  of  waves  and  fluid-mud  density  changes  were  obtained  simultaneously  with  tide  height 
and  suspended-sediment  data.  Fluid-mud  samples  were  taken  to  determine  mass  physical 
properties  of  the  bottom;  fathometer  profiles  were  run  to  determine  mudbank  morphology;  and 
current  profiles  were  taken  over  a tidal  cycle  to  determine  the  nearshore  residual  current. 


Previous  Studies 


Quantitative  field  data  on  process-form  interactions  in  muddy-coast  environments  are 
lacking  (McCave,  1972;  Drake,  1976).  The  concepts  for  sediment  suspension  and  transport  by 
waves  that  were  derived  for  low-turbidity  waters  and  sandy  bottoms,  although  abundantly 
present  in  the  literature,  cannot  always  be  applied  to  highly  turbid  waters  found  along  muddy 
coasts.  Much  of  our  knowledge  of  the  interaction  between  waves,  fluid  mud,  and  suspended 
sediment  stems  from  a few  field  studies,  mainly  by  Dutch  and  Russian  scientists,  beginning  in 
the  early  1960s. 

Delft  Hydraulics  Laboratory  (1962),  working  on  the  coast  of  British  Guiana,  has 
undertaken  perhaps  the  most  ambitious  field  study  yet  conducted.  An  important  conjecture 
with  respect  to  sediment  transport  was  that  fluid  mud  could  be  transported  en  masse  to  the 
northwest  by  wave  orbital  scour.  Large  mudbanks  spaced  every  30-60  km  along  the  coast  were 
believed  to  propagate  from  southeast  to  northwest  by  the  transport  of  fluid  mud.  Wave 
refraction  diagrams  for  idealized  mudbanks  indicated  that  maximum  convergence  of  wave 
energy  occurred  in  mudbank  "troughs."  Normally  a wave  convergence  would  be  expected  on 
mudbank  crests.  It  was  hypothesized  that  a shoreward  wave  drift  occurred  in  "troughs," 
whereas  a return  flow  was  present  over  mudbank  "crests."  Suspended-sediment  concentration 
was  found  to  increase  with  wave  orbital  velocity,  leading  to  a regional  distribution  of 
suspended  sediment  that  was  related  to  wave  activity. 

NEDECO  (1965),  working  in  a similar  environment  in  the  Gulf  of  Thailand,  found  that  in 
regions  of  fluid  mud,  waves  were  attenuated  rapidly  as  a result  of  bottom  friction.  It  was 
shown  that  sediment  was  brought  into  suspension  by  waves,  then  transported  great  distances  by 
coastal  currents. 

Zenkovitch  (1967),  summarizing  work  by  Russian  scientists  in  the  Gulf  of  Po  Hai, 
reported  that  during  wave  disturbances  bottom  turbidity  increased  and  a "creamy  mass"  was 
formed.  Nearshore  bottom  water  became  most  turbid  when  depths  were  "a  few  tens  of  cm" 
and  waves  impinged  directly  on  banks  of  fluid  mud.  Waves  and  currents  were  found  to  remove 
complete  layers  of  soft  mud,  presumably  during  storms,  which  later  would  be  redeposited. 

NEDECO  (1968)  conducted  an  extensive  hydraulic  investigation  along  the  Surinam  coast 
for  navigational  purposes.  Fluid  mud  was  found  to  damp  wave  motion  and  create  conditions 
favorable  for  further  sedimentation.  The  erosion  of  coastal  muds  was  generally  attributed  to 
waves,  but  the  large-scale  mud  transport,  estimated  at  2.5  x 10®mVyr,  was  believed  to  require 
a residual  current  farther  offshore. 

Nair  (1976),  elaborating  on  an  earlier  study  by  Moni  (1970),  documented  the  existence  of 
mudbanks  along  the  southwest  coast  of  India  that  were  nearly  identical  in  morphology  and 
sedimentology  to  those  of  northeast  South  America.  These  mudbanks,  although  formed  only 
during  monsoon  seasons,  served  as  a storehouse  for  littoral  sediments  and  provided  a buffer  to 
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wave  attack.  It  is  interesting  to  note  that  boundaries  of  dangerously  shoal  accumulations  of 
fluid  mud  in  southwest  India  are  usually  delineated  by  visually  noting  the  extent  of  the  calming 
effect  on  incident  waves. 

A mathematical  model  to  describe  the  effect  of  a fluid-mud  bottom  on  shallow-water 
waves  was  given  by  Gade  (1958)  to  explain  the  loss  in  wave  energy  over  the  "mud  hole,"  a 
localized  region  of  fluid  mud  off  the  central  Louisiana  coast  (Morgan  et  al.,  1953).  Wave 
energy  loss  into  a bottom  that  behaves  as  a viscous  fluid  was  found  to  substantially  reduce 
wave  height;  this  energy  loss  greatly  exceeded  the  dissipation  over  a rigid  sand  bottom.  In  a 
more  general  model  where  bottom  sediments  could  not  be  considered  as  a viscous  fluid,  Gade 
(1959)  found  that  energy  dissipation  took  place  as  the  bottom  deformed  elastically.  In  this 
case  sediment  structure  was  maintained  and  the  bottom  was  not  displaced  laterally.  Further 
work  by  Tubman  and  Suhayda  (1976)  has  shown  that  Mississippi  Delta  muds  respond  to  wave 
pressure  forces  elastically  and  thereby  dissipate  wave  energy  at  a rate  that  is  consistent  with 
field  measurements.  Vertical  bottom  movements  1-2  cm  high  were  found  to  dissipate  wave 
energy  at  rates  one  or  more  orders  of  magnitude  greater  than  those  from  percolation  or  normal 
frictonal  effects. 

In  a series  of  laboratory  experiments  Lhermitte  (1958;  1960)  found  that  wave  orbits 
could  penetrate  into  fluid  mud  to  substantial  depths.  With  mud  viscosities  less  than  20 
centipoises,  2-m-high  waves  of  8-sec  period  were  reported  to  cause  oscillations  down  to  20-m 
depth.  Later  laboratory  experiments  performed  by  Migniot  (1968)  and  Doyle  (1973)  verified 
that  waves  did  impart  orbital  motions  to  fluid  mud.  At  the  "mudline,"  orbits  were  elliptical, 
whereas  at  the  maximum  depth  of  movement  only  horizontal  motions  were  found.  Analytical 
results  presented  by  Bea  et  al.  (1975)  indicated  that  storm  waves  20  m high  and  14  sec  in 
period  could  cause  soil  movement  to  60  m below  the  mudline  and  a crest-to-trough  "mudline 
wave"  approximately  1.5  m high.  These  studies  lead  to  the  important  conclusion  that  net  drift 
resulting  from  non-closing  wave  orbits  can  cause  mud  transport  in  the  direction  of  wave  travel, 
provided  offshore  bottom  slop>es  are  gentle. 

Despite  the  great  number  of  investigations,  certain  problems  concerning  muddy-coast 
dynamics  remain  largely  unsolved;  1)  the  origin  of  mudbank  and  fluid-mud  deposits,  2)  the 
attenuation  of  waves  on  a fluid-mud  bottom,  3)  the  details  of  the  mechanism  of  fluid-mud 
movement,  and  4)  the  effect  of  fluid  mud  on  coastal  processes.  In  the  ensuing  paragraphs, 
parts  of  these  questions  will  be  answered. 


FIELD  AREA  AND  METHODS 


Field  experiments  were  conducted  between  20  September  and  8 October  1975  in  the 
nearshore  muddy  waters  of  the  central  Surinam  coast  to  obtain  simultaneous  data  on  waves, 
fluid  mud,  and  suspended  sediment  during  various  stages  of  the  tide.  An  instrument  check  and 
preliminary  measurements  were  made  in  a similar  environment  along  the  Louisiana  coast  in 
November  1974. 

The  country  of  Surinam,  situated  midway  between  the  Amazon  and  Orinoco  Rivers,  is 
representative  of  coastal  conditions  along  the  1600-km  coastline  of  northeastern  South 
America.  The  warm,  tropical  climate  is  controlled  by  the  northeast  trade  wind  system. 
Although  severe  storms  and  hurricanes  do  not  occur,  alternating  wet,  windy  and  dry,  calm 
periods  are  present  throughout  the  year.  September  and  October  are  typically  dry  months,  and 
relatively  calm  sea  conditions  prevail.  Beaches  are  composed  of  consolidated,  prograded 
mudflats  and  are  backed  by  mangrove  swamps.  Eight  mudbanks,  extending  5-10  km  alongshore 
and  2-5  km  offshore,  front  the  coast  every  30-60  km.  Accumulations  of  fluid  mud  1-2  m thick 
are  present  on  western  flanks  of  mudbanks  and  occasionally  on  mudbank  crests.  Pockets  of 
sand  and  shell  that  often  form  between  the  rhythmic  mudbanks  comprise  approximately  2 
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percent  of  the  total  volume  of  coastal  sediments  (Allersma,  1968). 

Data  were  taken  at  four  field  sites,  each  representing  a different  hydrologic  setting 
(Fig.  1):  subaqueous  mudbank  (site  1);  intertidal  mudbank  (site  2);  Surinam  River  entrance  (site 
3);  and  area  between  mudbanks  (site  4).  Depending  on  stage  of  the  tide,  water  depths  ranged 
from  1 to  3 m on  the  subaqueous  mudbank,  0 to  2 m on  the  intertidal  mudbank,  5 to  8 nti  in  the 
Surinam  River  entrance,  and  4 to  7 m between  mudbanks.  Fluid  mud,  herein  defined  as  a 
sediment-water  mixture  in  which  the  sediment  concentration  is  greater  than  10,000  mg/1 
(Krone,  1962),  was  present  on  the  subaqueous  and  intertidal  mudbank  sites  and  in  the  Surinam 
River  entrance.  Where  fluid  mud  is  present,  the  bottom  is  defined  as  the  surface  of  the  fluid 
mud  layer,  even  though  survey  instruments  may  penetrate  through  it  to  more  consolidated 
sediments  (Odd  and  Owen,  1972,  p.  183). 

Time  series  data  on  wave  characteristics  (wave  height,  H;  wave  period,  T)  were 
recorded  simultaneously  with  fluid-mud  density  changes  on  the  subaqueous  and  intertidal 
mudbank  sites  (1  and  2).  Wave  data  only  were  taken  in  the  area  between  mudbanks  where 
consolidated  mud  was  present.  These  data  were  taken  with  a wave/fluid-mud  pressure-sensing 
device  designed  and  constructed  at  Coastal  Studies  Institute  for  use  in  muddy  environments. 


Two  pressure  transducers,  fastened  securely  to  4-cm  diameter  pipe  that  was  driven  into 
the  bottom,  formed  the  sensing  unit  of  this  system.  Figure  2 shows  the  position  of  pressure 
transducers  relative  to  the  fluid-mud  bottom.  On  the  first  arrangement  (site  1)  transducer  A 
was  located  near  the  level  of  low  water  (LW)  and  transducer  B was  placed  8 cm  into  the  fluid 
mud;  on  the  second  arrangement  (site  2)  both  transducers  (A  and  B)  were  in  the  fluid  mud.  In 
most  cases,  surface-wave  oscillations  were  recorded  with  transducer  A and  density  changes  in 
the  fluid  mud  were  recorded  from  the  pressure  differential  between  transducers  A and  B. 

All  wave  and  fluid-mud  density  data  were  recorded  as  analogue  signals  on  a two-channel 
Gould  Brush  strip-chart  recorder.  Before  reaching  the  recorder,  signals  were  processed 
electronically  to  remove  high-frequency  noise  and  to  enhance  resolution.  Details  of  this 
instrument  and  its  use  are  given  in  Appendix  I. 

Water  samples  were  taken  at  each  field  site  in  a 1-liter  modified  Van  Dorn  bottle  at 
0.5-1-hour  intervals  for  suspended-sediment  concentration  and  particle  size  determination. 
Surface  samples  were  taken  at  the  two  field  sites  over  the  mudbank  (sites  1 and  2);  surface, 
mid-depth,  and  bottom  samples  were  taken  between  mudbanks  (site  4)  and  in  the  Surinam  River 
entrance  (site  3).  A total  of  82  samples  was  collected. 

A tide  record  was  obtained  for  the  duration  of  the  field  experiments  using  a battery- 
powered  pressure  transducer  water-level  gage  (fabricated  at  Coastal  Studies  Institute)  that 
was  fastened  to  a platform  in  the  Surinam  River  entrance  (Fig.  1).  Two  minutes  of  data  were 
recorded  continuously  on  a Rustrak  recorder  every  30  min.  A 60-sec  time  constant  filter 
removed  high-frequency  surface  waves. 

Current  profiles  were  taken  in  the  area  between  mudbanks  (site  4)  during  a tidal  cycle 
using  a Bendix  Q-15  meter  (Bendix  Environmental  Sciences  Division).  Logistics  problems 
precluded  obtaining  current  data  in  the  shoal  water  over  banks  of  fluid  mud. 

A variety  of  background  data  was  also  obtained;  salinity  and  temperature  hydrocasts 
were  taken  at  each  field  site  every  hour  using  a Beckmann  RS5-3  inductive  salinometer 
(Beckmann  Instruments,  Inc.);  eight  samples  of  fluid  mud  were  taken  over  various  parts  of  the 
mudbank  and  in  the  Surinam  River  entrance;  fathometer  profiles  were  run  to  de. ermine 
mudbank  morphology,  offshore  slope,  and  the  extent  of  fluid  mud;  and  an  extensix  aerial 
reconnaissance  was  undertaken  to  observe  coastal  landform  variability. 

All  time  series  data  (waves,  currents,  tide  height,  fluid-mud  density)  were  digitized  on  a 
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Calmagraphic  III  disc  digitizing  system  (Calma  Corporation).  A standard  stochastic  analysis 
was  performed  on  wave  records  using  a Biomedical  Time  Series  Program  (see  Appendix  II). 
Current  data  were  time  and  depth  corrected,  then  separated  into  tidal  and  nontidal 
components. 

Water  samples  were  filtered  through  Millipore  0.45-micron  filters  using  a pressure 
filtration  system  and  following  standard  procedures  outlined  by  Meade  et  al.  (1975).  Salinity 
was  determined  from  the  filtrate  of  each  sample  using  an  Elco  Lab  Line  (Lab  Line  Instruments) 
conductive  salinometer. 

Mass  physical  properties  (sediment  concentration,  bulk  density,  water  content,  voids 
ratio)  were  determined  in  the  laboratory  for  each  fluid-mud  sample;  viscosity  was  determined 
graphically  on  the  basis  of  salinity  and  sediment  concentration  (Delft  Hydraulics  Laboratory, 
1962).  A standard  pipette  analysis  (Folk,  1968)  was  performed  on  samples  selected  from  each 
hydrologic  environment.  Laboratory  methods  used  in  the  above  analyses  are  discussed  in  detail 
in  Appendix  III. 


RESULTS 


Coastal  Morphology  and  Hydrograohy 

Accumulations  of  fluid  mud  appear  as  partial  or  "false  bottom"  returns  on  fathometer 
records  (Fig.  3).  Profiles  taken  over  Warappabank  (Fig.  1)  indicate  that  the  bottom  slopes 
gently  seaward  at  0.00047  to  0.00073  (0°1'36"  to  0®2'30").  Coastal  mudbanks,  each  similar  in 
form  to  Warappabank,  were  visible  from  the  air  at  low  tide  along  the  Surinam  coast  and, 
together  with  the  trough  areas  between  mudbanks,  displayed  the  spectacular  pattern  of 
erosion-accretion  already  well  documented  in  the  literature.  Approximately  75  percent  of  the 
coastline  is  estimated  to  be  stable  or  accreting.  An  example  of  shoreline  progradation,  as 
evidenced  by  stage  of  mangrove  growth,  is  shown  in  Figure  4. 

Warm,  saline  water  is  present  along  most  of  the  Surinam  coast.  In  situ  temperature  and 
salinity  readings  showed  little  variation  over  tidal  cycles  or  between  field  sites.  Within  5 km 
of  the  shoreline,  in  September  and  October,  temperatures  ranged  from  28.7®  C to  3I.2°C  and 
salinity  ranged  from  34.34  °/oo  to  35.86  °/oo.  Discharge  from  the  Surinam  River  (400  m’/sec) 
turns  rapidly  to  the  west  and  does  not  affect  temperature,  salinity,  or  coastal  processes  at  the 
field  sites  to  the  east. 

Characteristics  of  the  tide  are  important  in  fluid-itiud  dynamics.  Water-level  changes, 
taken  continuously  during  the  field  experiments,  are  shown  in  Figure  5a.  Major  features  of  the 
tide  are  1)  a well-behaved,  semi-diurnal  periodicity  and  2)  a moderate  to  high  range  (3.2  m at 
spring  tide).  Because  of  the  tide  range  and  gentle  bottom  slope,  fluid  mud  may  be  exposed  to  a 
distance  offshore  of  more  than  3 km  during  each  low  water.  An  estimated  25  km*of  fluid  mud 
on  Warappabank  is  exposed  twice  daily;  upon  exposure,  pore  waters  are  expelled  and  initial 
stages  of  dewatering  begin  (Fig.  5b). 

Tidal  currents  measured  at  site  4 in  the  area  between  mudbanks  are  rotary,  opening 
clockwise,  with  maximum  near-surface  and  near-bottom  velocities  of  55  cm/sec  and  40 
cm/sec,  respectively  (Fig.  6a).  A residual  current,  setting  to  the  northwest,  is  superimposed  on 
the  cyclic  tidal  currents  (Fig.  6b).  This  residual  flow,  known  as  the  Guiana  Current,  originates 
as  an  extension  of  the  south  equatorial  current,  then  swings  to  the  northwest  and  parallels  the 
coastlines  of  French  Guiana,  Surinam,  and  British  Guiana  (Metcalf,  1968).  Maximum  speeds 
associated  with  this  current  are  reported  to  be  100-200  cm/sec  at  a distance  30-40  km 
offshore.  Maximum  residual  surface  velocity  at  site  4,  5 km  offshore,  was  20  cm/sec. 


6 


Figure  3.  Fathometer  profile  of  fluid-mud  bottom,  western  edge  of  Warappabank. 


Wave  Characteristics 

Examination  of  sixty-five  20-min  wave  records,  taken  over  five  sampling  days,  indicates 
that  waves  in  the  nearshore,  with  respect  to  wave  form,  are  of  two  types;  sinusoidal  and 
solitary-like.  In  water  5-8  m deep  between  mudbanks,  waves  are  generally  sinusoidal,  0.5-1. 0 
m in  height,  and  12-15  sec  in  period.  Assuming  that  linear  wave  theory  applies,  then 


L=^tanh(^)  (1) 

gives  wave  length,  which,  together  with  water  depth  in  the  ratio  h/L,  provides  a classification 
for  waves  based  on  depth  of  water  in  which  they  travel.  The  observed  h/L  ratios  ranged  from 
0.05  to  0.10,  indicating  that  these  are  intermediate  (transitional)  water  waves  (Coastal 
Engineering  Research  Center,  1973). 

A remarkable  change  in  wave  shape  takes  place  when  sinusoidal  waves  first  propagate 
as  shallow-water  waves  (h/L  < 0.05)  over  a fluid- mud  bottom.  Wave  records  taken  over 
mudbanks  between  HW  + 2 hours  and  LW  + 4 hours  show  that  wave  form  closely  approximates 
that  for  theoretical  solitary  waves.  Isolated  crests  are  separated  by  flat  troughs  lying  at  still- 
water  level.  Wave  period  is  12-16  sec  and,  in  water  1-3  m deep,  wave  height  is  0. 1-0.5  m. 
Crests  are  steep  and  symmetrical,  with  lateral  continuity,  estimated  from  aerial  recon- 
naissance to  be  more  than  1 km  (Fig.  7).  Examples  of  the  two  types  of  wave  forms  are  given  in 
Figures  Sa  and  b;  the  theoretical  solitary  wave  form  is  given  in  Figure  8c. 
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Figure  it.  Prograding  coastline  showing  three  levels  of  mangrove  growth,  eastern 
Surinam  coast. 


It  is  interesting  to  note  that  waves  become  more  like  true  solitary  waves  in  appearance 
as  sediment  concentration  in  water  increases.  Figure  9 shows  a sequence  of  wave  records, 
taken  on  the  intertidal  mudbank  (site  2),  that  illustrates  this  effect.  This  series  of  wave 
records,  spanning  3 hours,  shows  a marked  change  to  solitary-like  waves  with  flat  troughs 
identically  at  SWL  as  surface  suspended-sediment  concentrations  increase  from  2750  to  20,000 
mg/1.  On  intertidal  portions  of  mudbanks,  to  be  discussed  later,  the  concentration  of  sediment 
is  related  to  tide  height,  so  that  in  many  cases  solitary  wave  appearance  is  also  related  to  tide. 

At  LW,  when  banks  of  fluid  mud  are  exposed,  waves  spill  directly  onto  the  fluid  mud. 
During  a rising  tide  there  is  no  return  flow  associated  with  incoming  waves;  each  wave  simply 
advances  farther  toward  shore  and  the  water  motion  becomes  borelike.  On  several  occasions  it 
was  observed  that,  as  waves  impinged  on  a bank  of  fluid  mud,  wave  energy  was  transferred 
from  water  to  the  fluid  mud  itself  and  a solitary-like  wave  propagated  600  m through  fluid  mud 
to  the  shoreline. 

Field  observations  also  revealed  that,  except  during  LW,  when  waves  spilled  onto  banks 
of  fluid  mud,  little  wave  breaking  took  place.  Rather,  wave  height  was  continuously 
attenuated  as  waves  traveled  with  a solitary-like  wave  form  over  1-3  km.  Most  waves  do  not 
reach  the  shoreline.  A plot  of  wave  height  versus  water  depth  for  25  wave  records  taken  over 
the  mudbanks  indicates  that  over  fluid  mud  waves  decrease  in  height  linearly  with  decreasing 
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Figure  5b.  Photograph  illustrating  mud  exposure  at  low- 
water  spring  tide  in  the  vicinity  of  field  sites  1 and  2:  (A) 
muddy  coastal  water,  (B)  bank  of  exposed  fluid  mud,  and 
(C)  newly  accreted  wedge  of  fluid  mud  forming  a coastal 
mudflat  that  is  being  actively  colonized  by  mangroves. 
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Figure  5a.  Tide  record  in  the  Surinam  River  entrance  during  the  period  of  inves 
tigation.  This  gage  was  not  tied  in  to  a geodetic  or  leveling  net. 
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Figure  7.  Solitary-like  waves  from  an  altitude  of  250  m. 
Note  overall  crestal  continuity  (center)  and  disappearance 
of  waves  in  fluid  mud  (upper  right). 


Figure  8.  Sinusoidal  wave  form  in  area  between  mudbanks  (a),  solitary-like 
wave  form  over  subaqueous  mudbank  (b),  and  wave  form  as  given  by  solitary 
wave  theory  (c). 
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wave  form  with  increasing  sediment  concentration  and  decreasing  water  depth,  intertidal  mudbank 


water  depth  (Fig.  10).  An  interesting  feature  of  this  figure  is  that  the  slope  of  the  least 
squares  best  fit  line  corresponds  to  an  H/h  of  0.23;  96  percent  of  the  H/h  ratios  measured  were 
less  than  0.33.  The  H/h  ratio  at  which  waves  break,  from  solitary  wave  theory,  was  shown  by 
McCowan  (1894)  to  be  0.78.  Laboratory  and  field  measurements  over  rigid  bottoms  show  H/h 
ranges  between  0.73  and  1.03.  Thus  as  solitary-like  waves  advance  into  shoaling  water  over  a 
fluid-mud  bottom,  they  appear  to  have,  on  the  average,  an  H/h  ratio  considerably  less  than  the 
value  that  might  be  expected  from  solitary  wave  theory  or  from  measurements  over  sandy 
bottoms.  This  indicates  high  levels  of  wave-energy  dissipation. 

Additional  information  on  shallow-water  solitary-like  waves  can  be  obtained  from  the 
wave  spectrum.  Typical  spectra  of  waves  at  site  1 (subaqueous  mudbank)  are  given  in  Figure 
11.  Three  features  are  noteworthy  in  the  spectra.  First,  these  spectra  contain  a low- 
frequency  energy  peak  in  the  infragravity  range,  33-80  sec  (0.03-0.0125  cps).  Next,  there  is  a 
wide  distribution  of  wave  energy  over  all  frequencies;  on  many  of  the  spectra  examined  these 
secondary  energy  peaks  occur  at  harmonic  frequencies.  Third,  the  decay  rate  of  secondary 
peaks  is  roughly  exponential. 

The  effects  of  a fluid-mud  bottom  on  waves  in  shallow  water,  as  manifested  in  the  wave 
spectrum,  can  be  observed  in  Figure  12a-d.  This  figure  reveals  some  interesting  points.  A 
"characteristic"  energy  peak,  encompassing  a narrow  range  of  frequencies,  is  present  in  the 
spectrum  of  waves  from  the  area  between  mudbanks  (Fig.  12a).  Only  1 percent  of  the  wave 
energy  appears  to  be  contained  in  frequencies  higher  than  0.2  cps,  although  the  natural 
filtering  effect  of  the  pressure  transducer,  placed  4 m below  the  surface  in  this  instance, 
accounts  for  some  of  the  apparent  high-frequency  loss.  The  next  two  spectra  (Fig.l2b-c)  are 
of  waves  moving  across  identical  bottom  slopes  (3  = 0.0005),  each  in  1.3-m  water  depth,  but 
over  different  bottom  consistencies.  Waves  on  the  Louisiana  coast,  where  the  mud  is  more 
consolidated,  have  a sinusoidal  form  and  relatively  narrow  range  of  frequencies,  whereas  those 
on  the  Surinam  coast  are  solitary-like  and  have  an  apparently  wide  distribution  of  energy. 
Some  similarity  in  the  distribution  of  wave  energy  to  a wave  spectrum  generated  from  perfect 
solitary  waves  is  shown  in  Figure  12d.  The  difference  in  vertical  scales  should  be  noted. 

Although  the  peak  energy  density  associated  with  the  waves  between  mudbanks  is  more 
than  an  order  of  magnitude  greater  (approximately  25x)  than  for  waves  over  the  mudbanks, 
root  mean  square  height  is  only  1.5  times  greater,  whereas  total  wave  energy  is  3 times 
greater. 

Angle  of  wave  approach  at  the  shoreline,  as  determined  from  visual  estimates  in  the 
field  and  from  measurements  on  aerial  photographs,  ranges  from  3“to  S'in  shallow  water  over 
mudbanks  and  5®  to  10*  between  mudbanks. 


;nded  Sediment 


The  highest  overall  suspended-sediment  concentrations  and  the  greatest  variability  in 
concentration  occur  over  the  mudbanks,  whereas  the  lowest  and  least  variable  concentrations 
occur  between  mudbanks.  These  results,  plus  data  from  the  Surinam  River  entrance,  are  given 
in  summary  form  in  Table  1. 


Figure  13a-d  shows  details  of  suspended-sediment  change  with  stage  of  the  tide.  The 
major  features  of  concentration  change  over  the  subaqueous  mudbank  are  a rapid  increase  at 
mid-tide,  peak  at  low  water,  then  steady  decrease  to  the  following  high  tide.  A concentration 
change  of  more  than  two  orders  of  magnitude  takes  place  in  surface  water  during  each  tidal 
cycle.  On  the  intertidal  portion  of  the  mudbank,  concentration  changes  are  even  more 
pronounced.  Because  sampling  began  precisely  at  LW,  when  fluid  mud  was  exposed,  the  initial 
concentration  of  approximately  200,000  mg/1  shown  in  Figure  13b  is  that  of  fluid  mud  rather 
than  of  suspended  sediment.  As  water  level  rose,  these  suspended-sediment  concentrations 
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Figure  10.  Wave  height  versus  water  depth  for  twenty-five  20-min  wave 
records  taken  at  the  subaqueous  and  intertidal  mudbank  field  sites. 


decreased.  The  relationship  between  tide  height  and  concentration  for  all  samples  collected 
over  the  mudbank  is  given  in  Figure  14. 

In  contrast  to  the  sites  described  above,  no  trend  can  be  observed  with  stage  of  the  tide 
in  the  region  between  mudbanks  (Fig.  13d),  and  concentrations  remain  1-2  orders  of  magnitude 
lower  than  elsewhere.  Regardless  of  stage  of  tide,  concentrations  increase  with  depth. 

In  the  Surinam  River  entrance,  a region  also  influenced  by  fluid  mud,  sufficient  data 
were  not  available  to  determine  the  presence  of  a well-defined  trend.  As  before,  suspended- 
sediment  concentrations  increase  with  depth,  and  in  the  region  near  the  bottom,  fluid  mud 
persisted  at  all  stages  of  the  tide. 


Fluid-Mud  Dynamics 

Table  2 gives  characteristic  values  of  bulk  density,  water  content,  voids  ratio,  and 
viscosity  for  fluid  mud  in  several  locations  along  the  coast  and  one  location  in  the  Surinam 
River.  Bulk  density  is  always  highest  in  fluid  mud  that  is  undergoing  dewatering,  as,  for 
example,  on  the  intertidal  mudbanks,  which  are  exposed  1-3  hours  during  each  LW.  Minimum 
bulk  density  was  observed  in  fluid  mud  in  the  Surinam  River  channel,  where  strong  currents 
prevent  consolidation.  It  should  be  noted  that  considerable  variation  in  these  values  within  a 
given  fluid-mud  body  is  possible  since  densities  increase  near  the  base  of  the  fluid-mud  deposit. 
Water  content  and  bulk  density  are  inversely  related,  but  in  a strongly  nonlinear  fashion. 
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Figure  11.  Typical  wave  spectra  from  the  subaqueous  mudbank  field  site. 


Viscosity  increases  exponentially  with  bulk  density  and  spans  over  four  orders  of  magnitude 
within  the  range  of  fluid-mud  concentrations  observed  in  this  study. 

Typical  particle-size  distributions  for  coastal  and  river  muds  are  given  in  Figure  15. 
Mqdian  diameter  of  particles  in  unpeptized  form  is  0.5  to  1.0  micron,  although  it  is  likely  that 
particle  aggregates  many  times  larger  can  exist.  Less  than  1 percent  of  individual  particles  is 
coarser  than  clay  size,  and  the  brown,  fine  muds  along  the  Surinam  coast  are  remarkably 
uniform  in  size. 
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Figure  12a-d.  Contrasting  wave  spectra  for  mud  bottoms.  Intermediate-water 
waves,  consolidated  mud  bottom  (a),  shallow-water  waves,  consolidated  mud 
bottom  (b),  shallow-water  waves,  fluid-mud  bottom  (c),  and  solitary  wave 
spectrum  from  theory  (d).  Note  the  difference  in  vertical  scales. 


Observations  made  on  a bank  of  fluid  mud  at  LVV  indicate  several  interesting  features. 
Within  an  hour  after  fluid  mud  is  exposed  on  an  ebbing  tide,  a layer  (1-3  cm  thick)  of  very  clear 
water  appears  presumably  owing  to  the  initial  stages  of  dewatering.  This  water  gives  the  fluid 
mud  a glassy  appearance  when  observed  from  the  air  (Fig.  5b).  Under  shear  from  a northeast 
wind,  estimated  at  15  kn,  particles  from  the  surface  of  the  fluid  mud  are  picked  up  and  carried 
in  suspension  within  this  water  layer  as  it  begins  to  flow  Bottleneck  flowage  and  tensional 
cracks  that  have  vertical  walls  several  centimeters  thick  provide  evidence  that  the  entire  bank 
of  fluid  mud  moves  when  it  is  exposed  at  low  tide.  The  thixotropic  character  of  the  fluid  mud 
was  revealed  when  small  mass  movement  failures  were  gently  disturbed;  the  fluid  mud  lost  its 
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Table  1 


k 


Maximum  and  Minimum  Suspended-Sediment  Concentrations  at 
Four  Field  Sites,  Central  Surinam  Coast 


Sampling 

Concentration 

Range  in  Water 

Hydrologic  Setting 

Depth 

Maximum  Minimum 

Depth 

(10^  mg/1) 

(m) 

Subaqueous  Mudbank 

Surface 

37.49 

0.34 

1-3 

(Site  1) 

Intertidal  Mudbank 

Surface 

2222.85* 

26.57 

0-2 

(Site  2) 

Surinam  River 

Surface 

11.98 

0.49 

Entrance 

Mid-depth 

33.02 

0.64 

5-8 

(Site  3) 

Bottom 

303.62* 

240.02* 

Between  Mudbanks 

Surface 

0.27 

0.14 

(Site  4) 

Mid-depth 

Bottom 

0.82 

1.22 

0.16 

0.92 

4-7 

^"Suspensions"  with  concentrations  greater  than  100  x 10^  mg/1  are 
generally  referred  to  as  fluid  mud  (see  Krone,  1962).  Maximum  values  at 
site  2 correspond  to  samples  from  mudbank  surface  exposed  at  low  tide. 


internal  strength,  which  was  originally  sufficient  to  maintain  a vertical  wall,  and  the  tensional 
features  collapsed  instantly. 

Fluid  mud  responds  rapidly  to  the  passage  of  waves.  Figure  16  shows  examples  of 
density  changes  that  occur  on  the  subaqueous  mudbank  in  a 40-cm  interval  near  bottom  For 
reasons  discussed  in  Appendix  1,  two  scales  are  present  on  each  differential  density  record,  one 
for  long-term  density  changes  and  the  other  for  wave-by-wave  changes.  As  each  wave  passes, 
fluid  mud  is  suspended  and  a density  change  of  0.025  to  0.05  gm/cm^  takes  place.  Careful 
examination  of  these  records  reveals  the  important  fact  that  maximum  density  precedes  the 
wave  crest  by  0.5-1. 0 sec.  At  the  instant  of  crest  passage,  density  is  approximately  one-half 
its  peak  value,  and  1-2  sec  after  the  crest  passes,  near-bottom  density  reaches  its  lowest 
value.  The  magnitude  of  wave-by-wave  density  variation  decreases  after  1255  hours,  even 
though  the  surface  waves  remain  nearly  the  same  height.  It  should  be  noted  that  values 
obtained  from  this  figure  are  for  average  density  changes  between  pressure  transducers. 
Absolute  near-bottom  density  could  not  be  obtained,  but  is  estimated  to  range  from  1.16 
gm/cm^  at  transducer  B to  1.03  gm/cm’at  transducer  A;  a sharp  gradient  occurs  in  the  region 
10-30  cm  above  transducer  B. 

These  wave-by-wave  changes  are  superimposed  on  a longer  term  differential  density 
variation  that  is  related  to  tide.  Beginning  at  1115  hours,  between  time  of  mid-tide  and  LW, 
two  changes  take  place.  First,  a slight  increase  in  differential  density,  0.012  gm/cm^  occurs 
between  1115  and  1245  hours.  Next,  a rapid  decrease  of  0.040  grn/cm^  takes  place  between 
1245  and  1255  hours,  the  approximate  time  of  LW.  By  the  time  of  termination  of  the 
experiment,  at  1445  hours,  differential  density  was  reduced  an  additional  0.020  gm/cml 

Figure  17  shows  density  changes,  beginning  at  LW,  in  the  upper  0.5  m of  fluid  mud  on 
the  intertidal  mudbank.  Five  time  series  sections,  each  approximately  20  min  long,  show  that 
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Figure  13a-d.  Suspended-sediment  variations  at  four  field  sites.  Subaqueous  mudbanU  fa'; 
intertidal  mudbank  (b);  Surinam  River  entrance  (c);  area  between  rnudbanks  (d). 


bulk  density  decreased  from  approximately  1.16  gm/cm^  to  1.05  gm/cm^  over  a 3-hour  period 
as  water  level  rose.  An  important  feature  of  these  data  is  that,  although  water  did  not  begin 
to  cover  the  fluid  mud  at  the  field  site  until  1335  hours,  a measurable  change  in  average  fluid- 
mud  density  had  occurred  prior  to  this  time.  Considerable  density  variation  in  the  frequency 
-ange  : 5 min  is  evident;  however,  regression  analysis  shows  that  each  section  of  time  series 
has  a statistically  significant  slope. 


DISCUSSION 

Solitary  Waves  in  a Viscous  Fluid 

The  effect  of  a fluid-mud  bottom  on  waves  has  been  noted  for  more  than  two  centuries 
(Bristow,  1938).  Although  waves  in  shallow  water  are  often  considered  as  having  some  features 
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Figure  I't.  Suspended-sediment  concentration  versus  water  elevation  due  to 
stage  of  the  tide  for  all  samples  taken  over  the  mudbank. 


Table  2 

Physical  Properties  of  Fluid  Mud,  Surinam  Coast 


Location 

Bulk  Density 
(gm/cm^) 

Water  Content 
(%) 

Voids  Ratio 

Viscosity 

(cmVsec) 

Mudbank  exposed 
at  LW  (Site  2) 

1.172 

380 

10.1 

210 

Mudbank  covered 
by  2 m of  water 
(Site  1) 

Ll^tO 

490 

13.0 

75 

Mudbank-water 

interf  ace 

1.052 

2030 

53.8 

0.06 

Surinam  River 

entrance 

1.036 

4350 

115.3 

0.02 
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Figure  15.  Typical  particle-size  distributions,  central  Surinam  coast. 


of  a solitary  wave,  this  study  shows  for  the  first  time  that  the  form  of  these  waves  moving 
over  fluid  mud  is  remarkably  close  to  that  given  by  solitary-wave  theory. 

Taking  water  surface  elevation  in  the  time-space  domain  to  be 

[Vl^3 


n = H sech 


20 


DIFFERENTIAL  DENSITY 


WAVES 


WAVES 


I4A5  Hourt 


DIFFERENTIAL  DENSITY 


0.13 
1 0A8 
a 0.04 


0.13 

^ 0.08 

E 

■ 0.04 


0.13 

0.08 

xv 

£ 0.04 


0.13 
5 OXM 

x<s 

% 0.04 


Figure  16.  Variations  in  near-bottom  fluid-mud  density,  subaqueous  mudbank  (site  1).  Data 
were  taken  between  HW  + 4.5  hours  and  LW  + 2 hours. 
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Figure  17.  Time  series  of  bulk  density  changes  in  upper  0.5  m of  fluid  mud,  beginning  at  LW.  Values  above  each  time  series 
section  (A-E)  represent  surface  suspended-sediment  concentration. 


a solitary  wave  is  described  that  has  cis  its  major  features  a symmetrical  crest  entirely  over 
SWL  and  a trough  that  corresponds  exactly  to  the  level  of  still  water.  This  wave,  unlike  waves 
given  by  Airy,  Stokes,  and  Cnoidal  wave  theories,  is  specified  completely  by  wave  height  and 
water  depth. 

The  primary  difference  between  waves  shown  in  Figures  8b  and  9 and  waves  given  by 
equation  2 is  in  the  troughs.  An  average  trough  level  is,  in  many  cases  over  the  mudbank,  very 
close  to  SWL,  but  high-frequency,  low-amplitude  waves  are  present  between  crests  and  the 
troughs  are  seldom  completely  flat.  Dispersive  waves  in  trough  regions  make  it  difficult  to 
form  a true  solitary  wave  in  nature  (Coastal  Engineering  Research  Center,  1973,  p.  2.59). 

The  disappearance  of  these  dispersive  waves  with  increasing  sediment  concentration 
(Fig.  9)  is  an  excellent  example  of  the  filtering  effect  of  high  sediment  suspension  in  water  and 
offers  evidence  that  a fluid-mud  bottom  and  high  sediment  concentrations  in  the  water  column 
do  indeed  lead  to  wave  transformation  to  solitary-liKe  waves.  It  can  be  argued  that  waves  in 
the  areas  between  mudbanks,  upon  reaching  the  same  depth  of  water  as  solitary-like  waves 
over  mudbanks,  would  have  similar  characteristics.  Two  lines  of  evidence  suggest  that  this  is 
not  the  case.  First,  observations  made  from  a boat  traveling  parallel  to  bottom  contours 
revealed  a sharp  change  in  wave  form  when  it  crossed  onto  a region  of  fluid-mud  bottom.  It  is 
estimated  that  boundaries  of  fluid  mud  can  be  delineated  within  10  m simply  by  visual 
observation.  Second,  shallow-water  wave  spectra  from  Surinam  and  Louisiana  differ 
considerably,  even  though  waves  from  both  regions  travel  for  at  least  30  km  over  a mud 
bottom  with  identical  bottom  slope. 

A rigorous  evaluation  of  the  effect  of  fluid  mud  on  wave  characteristics  is  a formidable 
problem  that,  as  far  as  the  author  is  aware,  has  not  been  undertaken.  Results  from  this  study 
revealed  at  least  two  factors  that  complicate  such  an  evaluation.  First,  the  bottom  is  not  rigid 
but  is  highly  mobile  and,  relative  to  a fixed  consolidated  mud  bottom  beneath,  will  change 
F>osition  as  fluid  mud  is  suspended.  Second,  the  fluid-mud  bottom  undergoes  density  changes 
that  display  frequencies  from  those  of  waves  (~10-15  sec)  to  that  of  the  tide  (~12.4  hr). 
Viscosity  of  water  near  bottom  is  not  constant  inasmuch  as  fluid  mud  is  periodically  suspended 
into  the  water.  Even  lacking  the  density  and  viscosity  changes  attributed  to  waves,  the  near- 
bottom situation  is  complicated  because  a sharp  density  gradient  occurs  in  a bottom  layer 
estimated  to  be  10-30  cm  thick. 

Solitary  wave  theory  was  derived  under  the  simplification  of  an  inviscid  fluid  and  the 
boundary  condition  of  a rigid  bottom.  As  mentioned  above,  neither  of  these  assumptions  holds 
true  for  waves  over  mudbanks.  Equations  that  govern  wave  motion  in  general  neglect  viscous 
forces,  thus  giving  a solution  with  non-zero  velocity  at  the  bottom.  Oscillatory  waves  on  an 
impermeable  sandy  bottom  have  a no-slip  condition  imposed  at  the  boundary  so  that  the 
effects  of  viscosity  are  felt  only  in  a boundary  layer  in  the  immediate  vicinity  of  the  bottom 
(-1  cm  or  less).  Therefore,  viscous  forces  are  usually  negligible,  and  this  explains  the  success 
of  wave  theories  that  were  develop)ed  for  inviscid  fluids  (Madsen,  1976). 

One  effect  of  high  viscosity  and  a fluid-mud  boundary  such  as  exists  in  Surinam  may  be 
to  enlarge  the  boundary  layer.  Given  boundary  layer  thickness  under  a periodic  wave 

5 = (3) 

'‘(a)  ' 

where V is  kinematic  viscosity,  viscosities  in  the  range  of  1-10  stokes  for  15-sec  waves  will 
increase  the  boundary  layer  from  1 to  perhaps  10-30  cm.  In  fact,  the  entire  fluid-mud  deposit 
may  act  as  a boundary  layer.  NEDECO  (1965)  has  suggested  that  where  thick  accumulations 
of  fluid  mud  are  present  fluid  flow  in  the  near-bottom  boundary  layer  may  be  laminar  and 
mdy  have  Reynolds  numbers  of  about  250.  Under  laboratory  conditions,  estimated  boundary 
layer  thickness  beneath  solitary  waves  is  given  by  Ippen  et  al.  (1955)  as 
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where  the  wave  Reynolds  number, 
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(5) 


Again  taking  viscosities  in  the  range  of  1-10  stokes,  R^  <10*  and  6>10  cm  are  attainable. 
Low  Re  and  boundary-layer  growth  resulting  from  the  predominating  viscous  forces  are 
associated  with  high  rates  of  energy  dissipation. 

The  most  important  effect  of  increased  viscosity  in  the  water  column  is  wave  height 
attenuation.  Ordinarily,  waves  increase  in  height  with  decreasing  depth  up  to  the  point  of 
breaking.  In  the  general  case  for  linear  wave  theory,  wave  energy  flux,  P = ECn,  is  assumed 
to  remain  constant  from  deep  to  shallow  water  so  that 


P = (ECn)  = ECn  = CONSTANT 
o 


(6) 


where  the  subscript  "o"  represents  initial  conditions  and  the  energy  transmission  ratio,  n,  is 
defined  as  1/2  £n^l,  depending  on  water  depth. 


Now,  considering  solitary  waves  in  shallow  water  where  the  energy  transmission  ratio 


IS  unity 


In  this  case  E is  given  as 


C E = CD  . 
o o 


and  C,  given  to  a close  approximation  by  Boussinesq  (1872),  is 

C = i/ghlTTiT/h)  = *^+h)  . 

Substituting  equations  8 and  9 into  equation  7 yields 


H_  _o 

H ^ h 
o 


.h  A Uh/h’’^ 


(7) 


(8) 


(9) 


(10) 


This  is  shown  by  Ippen  and  Kulin  (1955)  to  be  nearly  equivalent  to 


H ‘'h  / 


4/3 


(11) 
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Thus,  on  the  basis  of  constant  energy  considerations,  for  solitary  waves  an  increase  in  wave 
height  would  be  expected  with  decreasing  water  depth.  Laboratory  studies  have  verified  the 
form  of  equation  1 1 for  predicting  wave  height,  although  wave  tank  experiments  show  that 
actual  heights  are  less  than  predicted  by  this  equation  (Ippen  and  Kulin,  1955).  However,  field 
data  from  the  present  study,  given  in  Figure  10,  indicate  that  over  fluid  mud  just  the  opposite 
occurs. 


That  the  discrepancy  between  actual  and  predicted  wave  heights  for  a given  set  of 
conditions  is  substantial  is  shown  in  Figure  18.  For  example,  according  to  solitary-wave 
theory,  a wave  100  cm  high  at  the  seaward  edge  of  fluid  mud  in  a water  depth  of  5 m would 
attain  a height  of  200  cm  after  traveling  to  a depth  of  only  3 m.  Shortly  after  that,  the  wave 
would  break.  In  reality  a wave  100  cm  high  follows  the  H/h  = 0.23  line  and  at  3-m  depth  is  60 
cm  high,  less  than  1/3  the  predicted  height.  This  wave  would  probably  not  break,  but  would 
continue  to  decrease  in  height  up  to  the  shoreline.  Whereas  Ippen  et  al.  (1955)  found  that 
damped  solitary  waves  maintain  their  original  volume  by  becoming  flatter  in  profile,  no  such 
change  in  wave  shape  was  noted  in  this  study. 


Wave  energy  loss  or  wave  damping,  as  manifested  in  height  attenuation,  can  be  caused 
by  many  factors:  bottom  percolation;  free  surface  dissipation;  internal  friction;  bottom  or 
boundary -layer  friction;  and  dissipation  into  a fluid  bottom.  Energy  losses  caused  by 
percolation  are  negligible  on  muddy  coasts  since  clays  are  highly  impermeable.  Attenuation 
at  the  free  surface  is  also  negligible,  although  Hunt  and  Brampton  (1972)  have  shown  that, 
under  certain  conditions  with  a contaminated  surface,  viscous  damping  may  occur. 


In  general  the  energy  loss  in  clear  water  resulting  from  internal  friction  is  quite  small, 
as  shown  by  Hough  (1896)  and  Lamb  (1932,  art.  349).  However,  when  viscosity  of  water  is 
increased  significantly  with  the  addition  of  sediment  particles,  to  say  1 stoke  or  more, 
internal  friction  may  become  important,  especially  for  high-frequency  waves.  Neglecting 
bottom  effects,  wave  height  can  be  shown  to  decay  exponentially  as 


H = H^ 


(12) 


where  vis  kinematic  viscosity  and  t is  time  of  travel  beginning  at  t = 0 when  wave  height  is 
H^  (Keulegan,  1950).  Taking  wave  period  = 15  sec,  viscosity  = 1 stoke,  and  t = 1000  sec,  the 
time  necessary  for  a solitary  wave  (H-0.5  m,  h~2  m)  to  travel  a distance  of  5 km  (equation  9), 
it  can  be  shown  that  H = 0.99  H^;  essentially  no  wave  energy  is  lost  in  this  instance  to 
internal  friction.  Since  H varies  as  the  fourth  power  of  wave  period,  the  height  attenuation 
for  short-period  waves  becomes  more  important.  For  instance,  given  a l-sec  wave  under  the 
same  conditions  as  above,  H = 0.04  H , and  the  wave  essentially  disappears  upon  reaching 
shore  after  traveling  5 km,  the  width  or a typical  mudbank.  In  a natural  situation  l-sec  waves 
will  travel  somewhat  slower  than  15-sec  waves  and  therefore,  in  the  hypothetical  situation 
above,  may  be  entirely  dissipated  before  reaching  the  shoreline.  Calculations  using  equation 
12,  based  on  typical  wave  data  gathered  from  regions  of  fluid  mud,  suggest  that  energy  loss  is 
non-negligible  over  the  width  of  a mudbank  for  waves  up  to  3-sec  periods.  This  is  important 
in  explaining  loss  of  high-frequency  waves  with  increasing  sediment  concentration,  as  shown 
in  Figure  9. 


Keulegan  (1948)  has  determined  that  the  theoretical  rate  of  loss  of  energy  owing  to 
viscous  shear  beneath  a solitary  wave  on  a smooth  horizontal  surface  is 
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where,  again,  v is  kinematic  viscosity.  Since  the  rate  of  energy  loss  must  be  manifested  in  the 
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Seaward  edge  of  fluid  mud 


Figure  18.  Plot  of  predicted  and  observed  wave  height  versus  water  depth  for 
solitary  waves, 

amplitude  attenuation,  the  dissipation  can  be  written 


^ ^ dE  ^ „ 

dt  ^ dH  dx  dt  ^ dH  dx 


where  C is  given  in  equation  9.  Equations  8 and  lit  may  be  combined  to  yield 


dE- 

dt 


H 


1/2  dH  . 
dx 


(15) 


Ippen  et  al.  (1955)  have  shown  that  equations  13  and  15  may  be  alternatively  written 


(16) 


where  is  initial  wave  height  and  H is  wave  height  after  the  wave  has  traveled  a distance  of 
X.  Laboratory  experiments  have  verified  this  mathematical  form  for  rate  of  attenuation 
(Daily  and  Stephen,  1953b;  Ippen  et  al.,  1955).  It  is  important  to  note  that  rate  of  energy  loss 
resulting  from  viscous  friction  is,  to  a first  approximation,  the  same  as  that  resulting  from 
bottom  friction.  In  each  case,  rate  of  energy  dissipation  per  unit  area  is  proportional  to  the 
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product  of  shear  stress  and  horizontal  particle  velocity  at  the  boundary  (Putnam  and  Johnson, 
1949) 


Df  = kt^u^ 


(17) 


where  k is  a friction  coefficient.  Numerical  equivalence  can  be  obtained  by  selecting 
appropriate  frictional  coefficients. 

To  assess  the  effect  of  viscosity  on  wave  attenuation  caused  by  viscous  friction,  the 
rate  of  attenuation  (equation  16)  was  plotted  for  viscosities  of  0.01,  0.05,  0.10,  and  1.0  stoke 
for  a constant  water  depth  beginning  at  an  H/h  of  0.78  (Fig.  19).  Taking  a viscosity  of,  say, 
0.05  stoke,  a wave  in  2 m of  water  will  undergo  a relative  height  attenuation  of  from  0.78  to 
0.45  after  traveling  5 km.  In  a water  depth  of  2 m this  corresponds  to  a height  decrease  of 
from  1.56  to  0.90  m.  For  viscosity  of  1.0  stoke,  the  corresponding  height  decrease  is  from 
1.56  to  0.44  m.  The  range  in  determined  H/h  ratios  for  the  two  field  sites  (1  and  2)  on 
Warappabank  have  been  plotted  in  Figure  19  at  the  appropriate  x/h  locations,  assuming  a 
mudbank  width  of  5 km  and  an  average  depth  of  2 m.  These  data  indicate  that  kinematic  vis- 
cosities of  0.1- 1.0  in  the  bottom  boundary  layer  are  sufficient  to  attenuate  waves  from  an 
H/h  = 0.78  to  the  H/h  ratios  observed  in  this  study.  If  initial  H/h  is  less  than  0.78,  a more 
realistic  condition,  then  a smaller  viscosity  is  required  to  achieve  the  observed  rate  of 
dissipation.  This  effect  is  clearly  significant  and  may  explain  in  part  the  steady  decrease  in 
wave  height. 

Although  equation  16  was  derived  for  a smooth  horizontal  bottom,  its  use  is  not 
restricted  to  laboratory  conditions.  The  nearshore  region  of  the  Surinam  coast  may  be 
considered  physically  and  hydraulically  smooth.  The  fact  that  bottom  slope  differs  from  zero 
and  therefore  water  depth  does  not  remain  constant  over  mudbanks  at  a given  stage  of  the 
tide  will  normally  cause  wave  shoaling.  Thus  additional  energy  must  be  dissipated  by  higher 
near-bottom  viscosities  in  the  boundary  layer  or  by  other  mechanisms. 

Gade  (1959)  has  developed  a model  to  explain  loss  of  wave  energy  resulting  from 
dissipation  into  a fluid  bottom.  The  motivation  for  this  model  was  the  observation  that  fluid 
mud  in  the  "mud  hole"  of  the  central  Louisiana  coast  significantly  altered  wave  heights 
(Morgan  et  al.,  1953;  Gade,  1958).  Basically,  the  model  required  a wave  motion  within  the 
bottom  fluid  mud  that  lagged  the  surface  water  wave  motion.  Tubman  and  Suhayda  (1976) 
have  shown  that  an  oscillating  bottom  can  dissipate  energy  at  a rate  at  least  an  order  of 
magnitude  greater  than  that  resulting  from  bottom  friction  alone.  As  will  be  discussed  in  the 
next  section,  it  is  not  known  whether  fluid  mud  over  mudbanks  oscillates  vertically  while 
maintaining  its  cohesive  structure;  however,  clouds  of  sediment  are  believed  to  be  suspended 
with  each  passing  wave,  and  it  is  hypothesized  that  large  but  unknown  amounts  of  energy  are 
lost  in  suspending  this  sediment. 


Fluid-Mud  Dynamics 

Susp>ension  and  deposition.  The  rapid  increase  in  suspended-sediment  concentration 
shown  in  Figure  13a  results  from  the  dispersion  of  bottom  fluid  mud  into  the  overlying  water 
column,  beginning  at  or  near  nid-tide.  Rata  plotted  in  Figure  14,  taken  on  three  different 
sampling  days,  show  that  as  water  depth  decreases  as  a result  of  falling  tide,  susp>ended- 
sediment  concentration  increases  exponentially.  The  coefficient  of  determination,  r*,  in- 
dicates that  60  percent  of  the  variation  in  susp>ended  sediment  can  be  attributed  to  changes  in 
tide  height.  The  unexplained  variability  may  be  due  to  height  of  incoming  waves.  Solitary 
wave  theory  predicts  that,  even  for  a constant  H/h,  higher  waves  have  greater  bottom 
velocities.  Thus  at  any  given  stage  of  the  tide,  higher  waves  may  suspend  more  sediment;  as 
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H increases,  so  does  h,  since  more  fluid  mud  is  suspended,  and  the  ratio  H/h  remains  nearly 
constant  at  0.23.  The  increase  in  water  turbidity  by  larger  waves  has  been  verified  qualita- 
tively from  coastal  reconnaissance. 

As  water  level  rises  over  mudbanks  exposed  at  LW  (Fig.  13b),  the  surface  layers  of 
fluid  mud  are  suspended  by  waves.  The  continuous  suspension,  combined  with  continuous 
increase  in  water  level,  results  in  a dilution  in  concentration.  By  termination  of  the  experi- 
ment shown  in  Figure  13b,  approximately  0.8  m of  fluid  mud  had  been  suspended  and  the 
bottom  was  effectively  lowered  by  this  amount. 

In  estuaries  high  suspensate  concentrations  may  also  be  attributed  to  the  presence  of 
fluid  mud  such  as  was  found  in  the  Surinam  River  entrance  (Fig.  13c).  Each  time  a mudbank 
migrates  past  the  river  mouth,  fluid  mud  gains  entrance  to  the  river  and  accumulates  in  the 
channel  as  lenses  1-2  m thick.  Strong  tidal  currents  sweep  mud  out  during  ebb  flows,  but  it 
returns  during  flood  and  concentrations  may  remain  high  for  several  years  during  the  passage 
of  a mudbank.  During  this  process  some  of  the  total  mud  volume  is  lost  to  tidal  overbank 
deposition.  Examination  of  water  quality  data  as  far  back  as  1961  in  various  reaches  of  the 
river  channel  (Hydraulics  Research  Division,  Ministry  of  Public  Works  and  Traffic, 
Paramaribo,  1970)  reveals  that  river  water  contains  average  suspended-sediment  concen- 
trations ranging  from  60  mg/1  to  22,000  mg/1  and  that  the  highest  of  these  values  occurred 
from  1966  to  1968,  a time  that  corresponded  with  the  last  passage  of  a mudbank  across  the 
river  mouth.  Thus  as  mudbanks  migrate  westward  each  river  entrance  or  harbor  is 
periodically  inundated  with  mud  that  originates  as  a coastal  deposit. 

The  lack  of  fluid  mud  in  areas  between  mudbanks  accounts  for  both  the  overall  low 
concentrations  and  the  lack  of  systematic  variations  in  concentration  (Fig.  13d).  These  low 
values,  15-80  mg/1  in  surface  waters  (Table  1),  are  actually  quite  typical  of  concentrations 
observed  on  other  muddy  coasts  throughout  the  world  (McCave,  1972).  The  high  concen- 
trations in  nearshore  water  over  mudbanks  are  unusual  and,  as  noted  by  Wells  and  Coleman 
(1977),  this  water  may  be  the  muddiest  in  the  world. 

The  suspension  of  fluid  mud  takes  place  in  two  modes,  as  indicated  in  Figures  16  and 
17:  short-term  wave-by-wave  suspension  and  longer  term  suspension  related  to  tide-induced 
water-level  fluctuations.  The  wave-by-wave  density  changes  in  Figure  16  are  caused  by  the 
suspension  of  fluid  mud  with  each  passing  wave.  A cloud  of  sediment,  with  concentration 
estimated  to  be  50,000  mg/1,  causes  these  density  changes,  as  shown  diagrammatically  in 
Figure  20.  At  time  t.,  8 cm  of  fluid  mud  covers  pressure  transducer  B.  A relatively  sharp 
density  gradient  extends  another  20  cm  upward,  where  density  decreases  from  1.16  gm/cm^  to 
1.026  gm/cm^,  the  density  of  seawater.  At  time  t2  , the  bed  is  instantly  suspended  and  a 
dense  cloud  of  muddy  water  extends  to  the  surface.  ^ These  sediment  clouds  were  not  visible 
during  the  experiment  because  of  the  ambient  high  turbidity.  This  suspension  causes  the 
average  density  between  pressure  sensors  A and  B to  decrease.  The  trend  toward  increasing 
differential  density,  beginning  at  t-,  results  from  the  gradual  settling  after  the  wave  crest 
has  passed.  Velocities  given  by  solitary  wave  theory  (Munk,  194^  indicate  that  water 
particles  first  move  at  a distance  of  x = lOh  from  the  crest,  but  that  particle  velocities  are 
only  10  percent  of  maximum  values  outside  the  crest  region  between  t2  and  t,.  This  may 
explain  the  suddenness  of  the  density  changes,  all  of  which  occur  within  a 2- 3-sec  period. 
Time  t.  represents  the  return  to  initial  bed  state;  near-bottom  density  increases  as  particles 
settle  during  the  time  between  passage  of  wave  crests.  Difficulties  in  explaining  this  high 
rate  of  settling  will  be  discussed  later. 

The  long-term  density  change  at  this  field  site  (Fig.  16)  is  also  quite  important  because 
it  shows  that,  simultaneous  with  the  suspension  of  clouds  of  sediment  by  waves,  a layer  of 
fluid  mud  is  slowly  stripped  from  the  bottom  near  time  of  low  tide.  This  is  the  layer  of  fluid 
mud  that  causes  the  high  suspended-sediment  concentrations  described  previously.  Between 
1245  and  1445  hours  (Fig.  16),  an  estimated  10  cm  of  fluid  mud  is  removed  from  the  bottom, 
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Figure  20a-b.  Wave-by-wave  near-bottom  density  changes  as  explained  by  sediment  suspension  and  deposition. 


much  of  it  during  the  10  min  from  1245  to  1255  hours.  The  substantial  decrease  in  wave-by- 
wave suspension  of  fluid  mud  at  1445  hours  is  caused  by  the  removal  of  this  fluid  mud  layer 
over  pressure  transducer  B.  Once  this  fluid  mud  is  removed,  pressure  transducer  B records 
very  nearly  the  same  density  changes  as  pressure  transducer  A,  thus  giving  a stable 
differential  density  record. 

Similarly,  the  dramatic  density  changes  on  the  intertidal  mudbank  (Fig.  17)  are  impor- 
tant in  understanding  fluid-mud  dynamics.  Between  the  hours  1210  and  1335,  before  water 
has  covered  the  mudbank,  both  long-  and  short-period  changes  in  fluid-mud  density  occur. 
These  density  changes  are  well  within  the  limits  of  precision  of  the  pressure-sensing  instru- 
ment (see  Appendix  I)  and  provide  evidence  that  movement  of  fluid  mud  takes  place  while  the 
mudbank  is  exposed  as  a tidal  flat.  The  propagation  of  solitary-like  waves  through  fluid  mud 
and  the  transport  of  clay  particles  in  clear  water  overlying  the  fluid  mud  are  other  indicators 
that  intertidal  banks  of  fluid  mud  are  not  stagnant  deposits. 

Another  important  feature  of  Figure  17  is  the  sharp  change  in  slope  between  time 
series  sections  D and  E,  at  a time  when  surface  waves  were  10  cm  high.  Aiiathurai  and  Krone 
(1975)  reported  that  at  bed  shear  stresses  just  above  the  critical  value  in  cohesive  sediments 
erosion  takes  place  particle  by  particle,  whereas  at  higher  stress  levels  the  bulk  shear 
strength  of  the  bed  is  exceeded  and  it  fails  totally,  resulting  in  instant  suspension.  The  rapid 
[ density  decrease  in  section  E may  correspond  to  such  a mass  erosion  phase  that  leads  to  rapid 

! and  complete  dispersion  of  fluid  mud  into  water.  Unfortunately,  no  data  are  available  after 

( mid-tide  and  further  details  of  this  process  cannot  be  determined.  It  is  known,  however,  that 

[ fluid  mud  is  redeposited  at  high  tide,  a conclusion  originally  based  on  the  clear  appearance  of 

I the  surface  water  and  later  verified  by  near-bottom  sampling  and  fathometer  profiling. 

Controversy  presently  exists  over  whether  fluid  mud  can  be  transported  as  a mass  by  a 
I creeping  motion  from  shear  exerted  by  overlying  water.  NEDECO  (1965)  and  Ippen  (1966) 

report  that  such  mass  movement  occurs  in  nature,  although  it  has  never  been  measured 
directly.  Parthenaides  (1971)  and  Krone  (1962;  1972)  feel  that  this  is  not  possible  and  that 
mud  must  first  be  suspended.  Evidence  from  visual  observation,  such  as  bottleneck  failures, 
indicates  that  mass  movement  without  suspension  takes  place  in  the  special  case  when  fluid 
mud  is  exposed  as  a LW  tidal  flat. 

The  related  question  of  whether  or  not  the  fluid-mud  bottom  oscillates  vertically  while 
maintaining  its  cohesive  structure  cannot  be  answered  from  this  study.  Documentation  of 
"mud  wave"  movements  in  soft  Mississippi  Delta  muds  (Suhayda  et  al.,  1976;  Tubman  and 
Suhayda,  1976)  suggests  that  this  is  likely  to  occur  in  Surinam;  however,  a "mud  wave" 
amplitude  of  1 5 cm  would  be  required  to  explain  the  wave-by-wave  density  changes  observed 
in  this  study.  A reasonable  assumption  is  that  some  vertical  movement  takes  place  near  the 
base  of  the  fluid  mud,  but  that  the  density  changes  can  be  caused  only  by  large  amounts  of 
mud  moving  as  a suspended  sediment  above  pressure  transducer  A (Fig.  20). 

Velocities  necessary  to  initiate  suspension  of  cohesive  sediment  have  generally  been 
determined  in  the  laboratory.  Considerable  discrepancy  exists  in  the  literature.  For 
example,  Zenkovitch  (1967)  reported  that  fine  silts  are  "stirred"  at  near-bottom  velocities  of 
7-12  cm/sec,  whereas  NEDECO  (1968)  reported  that  fluid  mud  is  not  suspended  until 
velocities  of  70  cm/sec  are  achieved.  Other  studies  show  that,  in  a single  series  of 
experiments,  a considerable  range  is  possible:  Allersma  et  al.  (1967),  20-100  cm/sec  for  Chao 
Phya  muds;  and  Delft  Hydraulics  Laboratory  (1962),  10-90  cm/sec  for  British  Guiana  muds. 
For  cohesive  sediments  in  general,  Drake  (1976)  suggested  that  average  velocities  of  10-30 
cm/sec  at  100  cm  above  the  bed  are  necessary  to  initiate  erosion.  Use  of  early  diagrams 
from,  Hjulstrom  (1935)  and  Sundborg  (1956),  as  reported  by  Parthenaides  (1971),  indicate 
values  of  18  cm/sec  for  erosion  of  unconsolidated  cohesive  sediments.  Although  no  current 
data  were  taken  over  banks  of  fluid  mud,  estimates  of  wave-induced  currents,  taken  from 
solitary  wave  theory,  reveal  the  large  magnitude  of  bottom  currents  even  for  small  waves 
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(Fig.  21).  On  the  basis  of  our  present  state  of  knowledge,  these  wave-induced  currents  from 
50  to  100  cm/sec  are  sufficient  to  suspend  fluid  mud. 

The  rapid  settling,  then  re-formation,  of  a fluid-mud  layer  during  a wave  period  is 
contradictory  to  one's  intuitive  expectation.  Generally  fluid  mud  forms  slowly  in  quiet  areas 
such  as  settling  and  turning  basins  or,  during  certain  stages  of  the  tide,  in  river  channels.  For 
individual  particles,  provided  that  concentrations  are  low,  settling  rate  is  given  by  the 
equation  for  Stokes  settling  velocity 


where  pis  dynamic  viscosity.  Median  particle  size  of  dispersed  clays  in  Surinam  is  approxi- 
mately 1 micron  or  less  (Fig.  15).  In  clear,  saline  water,  a 1-micron  clay  particle  will  settle 
at  1 X lO'"*  cm/sec  or  about  1 m in  12  days.  Clearly  this  settling  rate  is  many  orders  of  mag- 
nitude too  low  to  account  for  the  rapid  deposition  observed  in  this  study. 

In  water  where  salinity  exceeds  1 °/oo  , clay-sized  particles  aggregate  into  floes  due  to 
electrostatic  or  organic  bonding  (Whitehouse  and  Jeffrey,  1954;  Krone,  1962,  1972). 
Examination  of  Surinam  muds  verified  the  existence  of  such  floes;  Delft  Hydraulics 
Laboratory  (1962)  reported  that  floes  in  nearby  British  Guiana  muds  had  diameters  of  7-12 
microns.  The  settling  velocities  of  aggregates  rather  than  individual  particles  are  often  used 
to  explain  the  rapid  deposition  rates  in  muddy  regions.  Krone  (1972)  found  that  floes  in 
Savannah  Harbor  having  Stokes  settling  diameters  of  4-10  x 10^  microns  settled  from  0.45  to 
3.30  cm/sec.  According  to  Krone,  a 6 x 10^  micron  aggregate  contains  about  half  a million 
particles. 
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Figure  21.  Theoretical  bottom  particle  velocities  under 
solitary  wave  crests. 
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Settling  velocities  equal  to  those  observed  by  Krone  (1972)  in  Savannah  Harbor  are 
sufficient  to  explain  the  deposition  that  takes  place  after  sediment  is  suspended  by  the  tide. 
However,  these  rates  are  still  insufficient  to  explain  deposition  during  a wave  period.  An 
average  settling  rate  of  5-10  cm/sec  is  needed  for  the  rapid  mud  deposition  after  the  passage 
of  a wave  crest,  as  proposed  in  Figure  20.  This  is  the  settling  rate  of  quartz  spheres  1-2  mm 
in  diameter.  Quartz  sands  on  many  IJ.S.  East  Coast  beaches  are  smaller  than  this,  and  it  is 
difficult  to  visualize  mud  floes  of  this  equivalent  hydraulic  diameter. 

The  accepted  mode  of  settling  at  concentrations  greater  than  approximately  10,000 
mg/1  is  one  of  "hindered  settling"  whereby  interparticle  contact  is  maintained  and  water  is 
expelled  slowly  as  the  mass  settles  (Krone,  1962),  In  contrast  to  this  finding  in  the  literature, 
results  of  this  study  indicate  that  particles  may  settle  as  a unit,  a mass  that  fails  under  the 
shear  of  a passing  wave  but  rapidly  regains  its  structure  as  a fluid  mud  during  the  settling 
process,  then  forms  an  interface  that  quickly  collapses  to  give  an  overall  high  settling  rate. 
The  formation  of  such  mud-water  interfaces  during  settling  has  been  observed  by  several 
scientists,  most  notably  Pierce  and  Williams  (1966)  and  Owen  (1970).  Although  laboratory 
settling  rates  for  interfaces  are  ordinarily  no  more  rapid  than  for  "hindered  settling,"  field 
data  from  several  studies  suggest  that  settling  may  actually  be  faster  under  natural  con- 
ditions and  that  low-energy  environments  are  not  a prerequisite  for  rapid  settling.  For 
example,  Owen  (1969)  found  that  on  Thames  muds  in  situ  settling  velocities  were  an  order  of 
magnitude  higher  than  for  those  determined  under  quiescent  laboratory  conditions,  and 
McCave  (1971)  found  that  in  the  North  Sea  susp)ended-sediment  deposition  could  take  place 
regardless  of  the  value  of  a "wave  effectiveness  parameter"  for  sediment  scour. 
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Mud  transport.  The  important  finding  in  this  study  that  the  wave  profile  over  banks  of 
fluid  mud  is  very  close  to  that  given  by  solitary  wave  theory  has  been  utilized  to  assess  the 
role  of  waves  in  transporting  sediment.  Wave  height  and  water  depth,  together  with 
computed  wave  length  (equation  1)  for  each  wave  record  in  Figure  10,  show  that  waves  over 
mudbanks  are  clearly  within  the  application  bounds  for  solitary  wave  theory,  as  presented 
diagrammatically  by  Komar  (1976;  Fig.  22).  Although  only  wave  form  and  wave  period  have 
been  measured  in  the  field,  according  to  Munk  (1949)  a solitary-like  profile  alone  suggests 
application  of  this  theory.  This  is  perhaps  true  because,  unlike  most  wave  theories,  solitary 
wave  theory  was  developed  to  describe  what  had  been  observed  previously  in  nature.  Exten- 
sive laboratory  studies  by  Daily  and  Stephen  (1953a)  and  Ippen  and  Kulin  (1955)  have  shown 
that  if  wave  form  is  given  by  equation  2,  then  the  overall  behavior  of  the  waves  is  described 
by  solitary  wave  theory. 

Assuming  that  waves  over  mudbanks  can  best  be  described  by  solitary  wave  theory, 
then  the  potential  for  longshore  sediment  transport  is  high  for  three  reasons.  First,  hori- 
zontal and  vertical  water  particle  velocities  under  a solitary  wave,  given,  respectively,  by 


- r'M  l+cos(My/h)  cosh(Mx/h) 
[cos(My/hT+cosh(Mx/h^l  ^ 


(19) 


sin(My/h)  sinh  (Mx/h) 

~ [cos  (My/h)+cosh(Mx/h)  p 


(20) 


show  that  there  is  no  reversing  or  oscillatory  flow,  so  that  water  moves  only  in  the  direction 
of  wave  travel.  Second,  deepwater  wave  statistics  compiled  by  NEDECO  (1968)  from  40 
years  of  observations  indicate  that  93  percent  of  sea  and  swell  waves  approach  the  coast 
between  N30®E  and  east;  because  of  this  direction  of  wave  travel,  there  is  a strong  year- 
around  longshore  component.  Third,  as  given  in  Table  1,  suspended-sediment  concentration  is 
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Figure  22.  Generalized  wave  validity  diagram  showing  regions  of  applicability  of  several  wave 
theories  (from  Komar,  1976).  Note  region  of  waves  over  rnudbanks  in  Surinam. 


high  in  water  over  rnudbanks  where  solitary-like  waves  are  present.  Temporarily  ignoring 
other  coastal  processes,  the  net  outcome  will  be  a sediment-laden  water  that  moves  only 
shoreward  with  a strong  westward  component. 

In  previous  investigations  along  the  Surinam  coast  waves  were  considered  mainly  as  an 
important  eroding  and  suspending  agent  (NEDECO,  1968).  Sediment  transport  over  long 
distances,  such  as  between  the  Amazon  and  Orinoco  Rivers,  was  attributed  to  the  swiftly 
flowing  Guiana  current.  Estimates  of  sediment  transport  along  the  northeastern  coast  of 
South  America  are  given  variously  in  tons  and  cubic  meters  and  by  month  and  season.  These 
estimates,  taken  from  Delft  Hydraulics  Laboratory  (1962),  Gibbs  (1967),  Allersma  (1968), 
NEDECO  (1968),  and  Eisma  and  van  der  Marcel  (1971),  can  be  summarized  as  follows.  Ap- 
proximately 20-U0  percent  of  the  total  load  from  the  Amazon  River  is  transported  northwest 
along  the  coast.  This  amounts  to  150  x 10*  m’/yf  of  "through  transport"  in  the  form  of 
suspended  sediment  and  another  100  x 10*  mVyr  that  moves  in  the  form  of  mudbank 
propagation.  Thus  a total  of  250  x 10*  mVyr  moves  in  a band  that  extends  from  the  shoreline 
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to  30  km  offshore.  These  estimates  were  based  on  volume  and  rate  of  mudbank  migration, 
average  residual  current  speed,  and  suspended-sediment  concentration  on  the  continental 
shelf.  In  the  following  paragraphs,  volume  transport  that  can  be  explained  by  currents 
associated  with  solitary  waves  is  examined. 


Consider  mass  transport  to  be  given  by 


T 


M 


= cAv 


(21) 


where  c and  v are  time  and  depth  averages.  If  bulk  density  of  the  transported  mass  after 
deposition  is  known  or  assumed,  then  volume  transport  is  found  by  dividing  by  sediment 
concentration  at  that  density 


T = 

V 


(22) 


Next,  consider  transport  to  be  through  a cross-sectional  area  determined  from  the 
width  of  the  mudbank  and  average  water  depth  (Fig.  23).  Since  the  longshore  component  of 
transport  depends  on  angle  of  wave  approach,  equation  22  becomes 


T = cAv  sin  a . 

V p ~ • 

c 

Again  assuming  that  solitary  wave  theory  provides  a reasonable  description  of  waves  over 
mudbanks,  the  total  volume  per  unit  crest  length  above  SWL  can  be  obtained  as 


Figure  23.  Idealized  mudbank,  coast  of  Surinam.  Dashed  lines  represent 
longshore  components  of  wave  drift. 
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Q=  f r\dx  = 2 f 
-00  5 

where  the  substitution  of  equation  2 for  q yields 


Q = 2 / H sech' 

0 

Munk  (1949)  integrated  equation  25  to  obtain,  to  a good  approximation, 


Q = .h=  Vi  » . 


(25) 


(26) 


Since  a volume  transport  equal  to  Q takes  place  during  one  wave  period,  T,  the  mean  trans- 
port per  unit  time  equals  Q/T  and  volume  transport  velocity,  averaged  from  the  surface  to 
the  bottom,  yields 


V Q_ 

hT  ^ f 


(27) 


Equation  27  has  units  of  LT* ‘since  Q is  taken  to  be  volume  p)er  unit  length  of  wave  crest. 
Substitution  of  equation  27  into  equation  23  for  velocity  allows  a more  rapid  determination  of 
sediment  transport  since  particle  velocities  under  a solitary  wave  must  be  integrated 
numerically  for  each  wave  of  different  H/h  ratio.  Comparable  results  are  obtained  using 
either  particle  velocities  or  average  volume  transport  rate  of  water. 

Table  3 gives  values  of  longshore  transport  per  year  for  solitary  waves  of  15-sec  period 
arriving  at  the  shoreline  at  a 10°  angle.  Suspended -sediment  concentration  was  taken  to  be 
2000  mg/1,  averaged  over  a 5-km-wide  mudbank.  Transport  was  assumed  to  take  place  3 
hours  each  tidal  cycle,  1.5  hours  before  LW  and  1.5  hours  after  LW.  This  corresponds  to  the 
length  of  time  that  the  most  solitary-like  waves  were  observed  in  the  field. 

The  determined  volume  transport  increases  with  increasing  water  depth  and  wave 
height.  Figure  24  shows  this  rate  of  volume  transport  increase  assuming  constant  H/h  of  0.23, 
that  observed  in  the  field  (Fig.  10),  and  0.78,  the  limit  of  solitary  wave  theory.  A change  of 
approximately  two  orders  of  magnitude  takes  place  as  wave  height  increases  from  10  to  100 
cm,  maintaining  a constant  H/h  ratio.  The  lower  transport  rates  for  higher  H/h  ratios  (0.78 
versus  0.23)  are  due  to  the  fact  that,  for  given  wave  height,  water  depth  is  less  for  H/h  = 0.78 
than  for  H/h  = 0.23,  and  volume  of  water  in  a solitary  wave  is  related  to  the  square  of  water 
depth  (equation  26). 


Volume  transport  will  increase  linearly  with  increases  in  suspended-sediment  con- 
centration, cross-sectional  area,  and  portion  of  the  tidal  cycle  for  which  waves  are  considered 
solitary.  Higher  angles  of  wave  approach  will  also  lead  to  higher  transport  rates;  however, 
only  low  angles  of  approach  are  present  at  the  shoreline  in  Surinam.  The  variable  that  is 
perhaps  most  likely  to  affect  transport  is  suspended- sediment  concentration.  As  previously 
discussed,  the  vertical  distribution  of  suspended  sediment  is  poorly  known.  If  high  concen- 
trations, say  2000  mg/1  to  20,000  mg/1,  travel  near  bottom,  as  has  been  found  for  estuarine 
fluid  muds  (Inglis  and  Allen,  1957),  then  the  average  concentration  may  be  considerably  higher 
than  the  2000  mg/1  taken  here.  As  discussed  in  the  last  section,  clouds  of  sediment  with 
concentrations  of  50,000  mg/1  are  suspended  by  each  wave.  Much  of  this  will  be  transported. 
Further,  the  transport  calculations  are  based  on  wave  and  sediment  data  taken  during 
September  and  October,  the  two  calmest  months  of  the  year  in  terms  of  sea  conditions.  Data 
taken  during  rough-season  months  may  increase  these  values  significantly. 
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Figure  2k.  Volume  transport  of  sediment  by  solitary  waves 
for  H/h  = 0.23  (observed)  and  H/h  = 0.78  (theory). 


Table  3 

Longshore  Volume  Transport  of  Mud  (10'’m7yr)  by  Solitary  Waves 
As  a Function  of  Wave  Height  and  Water  Depth 


h(cm) 

H (cm) 

30 

40 

50 

60 

70 

70 

2.70 

3.12 

3.48 

» 

» 

90 

3.92 

4.52 

5.10 

5.58 

6.02 

110 

5.30 

6.12 

6.84 

7.56 

8.16 

130 

6.84 

7.58 

8.78 

9.66 

10.46 

*H/h  exceeds  0.78. 


NOTE;  Values  were  determined  for  T = 15  sec,  a = 10“  , 
c = 2000  mg/1,  W = 5 km. 
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Typical  volume  transport  rates  for  U.S.  East  and  West  Coast  sand  beaches  range  from  1 
to  3 X lO^mVyr  (Johnson,  1956;  1957)  and  a transport  rate  of  10®mVyr,  found  on  some  high- 
energy  California  beaches,  is  considered  unusually  high.  The  rates  of  3-70  x 10®  m’ calculated 
above  for  mudbanks  (Table  3;  Fig.  24),  even  as  order  of  magnitude  estimates,  are  perhaps  100 
times  higher  than  on  sandy  coasts.  If  the  volume  transport  estimates  by  Allersma  (1968)  and 
NEDECO  (1968)  are  applied  to  a 5-km  instead  of  30-km  width  of  mudbank,  then  the  volume 
transport  over  this  region  of  fluid  mud  would  be  approximately  42  x 10®  mVyr.  Thus  wave- 
induced  currents  can  readily  explain  these  high  rates  of  transport  without  the  aid  of  a 
superimposed  current.  Zenkovitch  (1967)  has  found  volume  transport  in  the  Gulf  of  Po  Hai  to 
be  10®  mVyr,  a fact  that  suggests  that  other  mud  coasts  have  similarly  high  rates. 

Although  volume  transport  may  be  many  times  higher  in  muddy  than  in  sandy  environ- 
ments, owing  to  the  high  water  content  of  fluid  mud  (Table  2),  mass  transport  may  actually  be 
less.  The  values  given  in  Table  3 and  Figure  24  were  determined  from  a transported  sediment 
concentration  of  250  kg/m’,  that  of  fluid  mud  (bulk  density  1.18  gm/cm’).  In  studies  on  sand 
and  coral  beaches,  volume  transport  is  calculated  from  an  immersed  weight  transport  rate 
that  is  determined  empirically  from  wave  characteristics  (Bagnold,  1963;  Komar,  1976,  p. 

206).  This  method  takes  into  account  the  density  of  sediment,  density  of  water,  and  a 
correction  factor  for  pore  space  in  beach  sand  so  that  mass  and  volume  are  related,  as  in  the 
present  study,  essentially  by  bulk  density. 

This  technique  cannot  be  applied  to  muddy  coasts  for  two  reasons.  First,  sufficient 
field  data  are  not  available  for  "calibrating"  the  relationship  between  wave  characteristics 
and  immersed  sediment  transport  rate  when  the  sediment  is  carried  in  suspension  as  clay 
aggregates.  Second,  as  previously  discussed,  waves  over  mudbanks  are  not  sinusoidal  and  the 
application  of  this  technique,  based  on  energy  considerations  from  linear  wave  theory,  may 
not  be  valid.  An  important  point  with  respect  to  volume  and  mass  is  that  volume  transport  on 
muddy  coasts  may  actually  be  more  important  than  mass  transport.  Since  fluid  mud  occurs 
regularly  as  large  banks  that  migrate  as  a body  with  delineable  boundaries,  it  is  volume  and 
not  mass  of  the  body  that  affects  coastal  processes. 

Several  simplifications  that  require  qualitative  evaluation  have  been  introduced  into 
the  transport  calculations.  The  longshore  component  of  net  drift  has  been  determined  under 
the  condition  that  no  coastal  boundary  exists.  This  is  clearly  unreasonable.  The  shoreward- 
directed  net  drift  component  requires  a compensating  return  flow  to  satisfy  continuity 
considerations.  However,  if  a uniform  return  flow  is  taken  perpendicular  to  the  coastline, 
which  is  common  procedure  for  solitary  wave  applications  (Munk,  1949),  then  longshore 
transport  is  unaffected.  An  alternate  method  for  return  flow  is  through  the  generation  of  rip 
currents.  There  is  some  field  and  laboratory  evidence,  cited  by  Munk,  that  in  the  case  of 
solitary  waves  a slight  but  continuous  seaward  movement  of  bottom  water  does  take  place, 
thus  suggesting  that  rip  currents  alone  do  not  fully  balance  the  shoreward  transport  of  water. 

Longshore  circulation  patterns,  in  which  rip  currents  may  be  a part,  have  also  been 
ignored.  Because  of  the  irregular  nearshore  topography  resulting  from  the  presence  of 
mudbanks,  waves  are  refracted,  giving  rise  to  localized  convergence  and  divergence  of  wave 
energy.  According  to  Allersma  (1968)  and  NEDECO  (1968),  convergence  of  wave  energy  is  in 
the  area  between  mudbanks,  whereas  divergence  is  over  mudbanks.  Based  on  this,  each 
concludes  that  shoreward  transport  of  water  is  between  mudbanks  and  seaward  transport  of 
water  is  over  mudbanks. 

Some  doubt  may  be  attached  to  these  conclusions.  The  normal  wave  refraction  pattern 
is  for  wave  convergence  over  offshore  shoals  and  shoreface-attached  ridges  (Goldsmith, 

1976),  even  if  their  trend  is  oblique  to  the  shoreline  such  as  along  the  Surinam  coast.  The 
refraction  diagrams  presented  by  Allersma  (1968)  and  NEDECO  (1968)  are  considerably  i 

iversimplified  and  are  constructed  for  an  idealized  mudbank;  different  results  can  be  ob-  j 

tained  by  selecting  different  wave  ray  spacings  and  angles  of  approach.  j 
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The  impwrtance  of  the  areas  of  convergence  is  that  they  give  rise  to  higher  waves. 
Longshore  currents  are  then  generated  from  the  longshore  gradient  in  wave  energy.  This,  in 
turn,  may  set  up  a circulation  cell  so  that  return  flow  takes  place  in  areas  of  wave  diver- 
gence. Even  though  waves  are  higher  between  mudbanks,  where  they  expend  their  energy  as 
plunging  breakers,  it  is  possible  that  a greater  wave  set-up  occurs  over  the  mudbanks,  owing 
simply  to  the  translatory  nature  of  solitary-like  waves.  If  this  is  the  case,  then  an  onshore- 
offshore  flow  pattern  opposite  to  that  proposed  by  Allersma  (1968)  and  NEDECO  (1968)  would 
follow.  Recent  work  by  Sonu  (1972)  on  the  west  coast  of  Florida  documents  a shoreward 
current  over  shoals,  where  waves  break  by  spilling,  and  a return  flow  in  troughs,  where  waves 
break  by  plunging. 

In  connection  with  nearshore  circulation,  currents  measured  5 km  offshore  at  field  site 
4 indicate  that  a substantial  residual  current  is  present.  Closer  inshore,  the  effects  of  the 
Guiana  current  are  unknown  but  are  believed  to  be  considerably  less,  especially  over 
mudbanks.  As  pointed  out  by  Komar  (1976,  p.  203),  the  belief  that  tidal  currents  and  ocean 
currents  approach  close  to  shore  and  were  primarily  responsible  for  longshore  transport  has 
been  replaced  by  the  concept  of  wave-induced  longshore  currents  as  being  the  chief  cause  of 
sediment  movement. 

Data  from  NEDECO  (1968)  indicates  that  at  a distance  of  20  km  offshore,  residual 
surface  currents  range  from  approximately  20  to  40  cm/sec,  depending  on  longitude  and 
season.  This  is  a distance  offshore  where  sediment  concentrations  are  low.  In  this  study  it  is 
believed  that  most  of  the  sediment  transport  takes  place  where  thick  accumulations  of  fluid 
mud  are  present.  The  complete  suspension  of  fluid  mud  to  a distance  offshore  of  perhaps  5 
km,  then  transport  by  waves,  may  carry  far  more  sediment  than  the  continuous  residual 
current  in  low-concentration  areas. 

Studies  by  Eisma  (1967)  and  Gibbs  (1975)  have  shown  that  mud  derived  from  the 
Amazon  River  is  trapped  in  the  nearshore  region  by  net  shoreward-flowing  bottom  currents. 
Once  in  this  coastal  region,  the  major  transport  may  be  due  to  wave-associated  currents  and 
circulation.  The  wave/fluid-mud  interaction  tends  to  be  self-perpetuating;  once  waves 
become  solitary-like,  they  transport  the  sediment  both  shoreward  and  alongshore,  so  that 
mudbanks  may  move  by  the  transport  of  fluid  mud. 


Implications  to  Other  Coasts 

"Mudflats"  and  "mangrove  flats"  represent  23  percent  of  the  total  coastline  in  North 
and  South  America  (Dolan  et  al.,  1972)  and  perhaps  a similar  percentage  worldwide.  At  least 
five  muddy  coasts  that  are  open  and  thus  exposed  to  incoming  waves  have  accumulations  of 
fluid  mud:  central  Louisiana  coast  (Morgan  et  al.,  1953);  northeastern  coast  of  South  America 
(Delft  Hydraulics  Laboratory,  1962;  NEDECO,  1968);  Gulf  of  Thailand  (NEDECO,  1965);  Gulf 
of  Po  Hai  (Zenkovitch,  1967);  and  southwestern  coast  of  India  (Moni,  1970;  Nair,  1976).  In 
addition,  many  estuaries  such  as  San  Francisco  Bay  (Krone,  1962)  and  Delaware  Bay 
(Simmons,  1966)  and  river  entrances  such  as  the  Mississippi  (Bates,  1953;  Fisk  et  al.,  1954; 
McClelland,  1967)  and  Severn  (Kirby  and  Parker,  1973)  are  partially  filled  with  tenacious  fluid 
mud.  The  new  information  gained  in  this  study  may  be  important  to  understanding  suspension, 
deposition,  and  transport  dynamics  in  these  and  other  coastal  regions. 

On  each  of  the  muddy  coasts  above  where  fluid  mud  is  present,  two  features  are 
common:  high  suspended-sediment  concentration  and  extraordinary  change  in  wave  form.  A 
major,  implication  of  this  study,  then,  is  that  if  waves  moving  over  a fluid-mud  boundary  have 
characteristics  of  solitary  waves,  then  sediment  transport  by  waves  alone  is  possible.  In  fact, 
inasmuch  as  solitary-like  waves  have  been  observed  where  suspensate  concentrations  are 
highest,  namely  over  fluid-mud  deposits,  then  in  each  of  these  areas  transport  of  sediment  by 
waves  may  exceed  transport  by  stronger  offshore  currents  where  sediment  concentrations  are 
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lower.  Thus  an  appropriate  explanation  for  sediment  transport  in  muddy-coast  environments 
may  require  a combination  of  wave-associated  currents  in  the  nearshore  fluid-mud  region  and 
a residual  oceanic  current  farther  offshore.  Furthermore,  high  volume  transport  rates,  say 
greater  than  10  x 10®  m^,  are  possible  without  breaking  waves  and  utilization  of  the  concepts 
of  radiation  stress  and  a nearshore  circulation  cell,  as  normally  applied  to  sandy  coasts 
(Komar,  1976,  p.  170),  even  though  such  patterns  of  circulation  may  occur. 

Before  the  importance  of  mud  transport  by  waves  can  be  further  evaluated,  detailed 
field  measurements  must  be  undertaken.  At  present  less  than  10  percent  of  the  continental 
shelf  areas  have  been  studied  for  suspended  sediment  (Drake,  1976)  and  detailed  wave  and 
current  data  are  available  for  an  even  smaller  percentage. 


CONCLUSIONS 


Analysis  of  field  data  has  shown  that  shallow-water  waves,  fluid  mud,  and  suspended 
sediment  interact  in  the  nearshore  region  of  the  Surinam  mud  coast.  These  interactions  are 
not  always  straightforward  and  may  be  contradictory  to  findings  in  sandy-coast  environments. 

The  presence  of  a fluid-mud  bottom  affects  incoming  waves  by  changing  their  form 
from  sinusoidal  to  solitary-like  arid  by  preventing  wave  breaking  except  for  occasional 
spilling.  Measurements  from  two  field  sites  indicate  that  as  solitary-like  waves  12-16  sec  in 
period  and  10-50  cm  high  propagate  shoreward,  their  wave  height  to  water  depth  ratio,  H/h, 
remains  nearly  constant  at  0.23,  regardless  of  water  depth.  At  higher  suspended-sediment 
concentrations,  troughs  become  flatter  and  the  solitary-like  apjjearance  becomes  more 
pronounced.  When  examined  by  Fast  Fourier  techniques,  periodic  solitary-like  waves  with 
flat  troughs  and  a sech^crest  profile  give  rise  to  a distribution  of  wave  energy  that 
characteristically  contains  harmonics. 

It  is  concluded  that  the  absence  of  high-frequency  waves  (T  < 3 sec)  in  solitary  wave 
troughs  may  be  due  to  internal  friction  from  the  increased  viscosity  of  water  by  the  addition 
of  sediment  particles.  The  steady  decrease  in  wave  height  without  breaking,  up  to  the  shore- 
line, is  explained  in  part  by  energy  dissipation  within  the  bottom  boundary  layer.  Substantial 
but  unknown  amounts  of  wave  energy  may  be  lost  by  bottom  oscillations  and  sediment  suspen- 
sion. An  important  conclusion  concerning  wave  attenuation  is  that  muddy  coasts  tend  to 
protect  their  shorelines,  particularly  the  fluid  mud  portions,  since  wave  energy  is  often 
dissipated  entirely  before  waves  reach  the  shoreline. 

As  solitary-like  waves  move  across  a mud  bottom,  the  soft  fluid  mud  is  suspended 
rapidly.  On  a wave-by-wave  basis,  near-bottom  measurements  indicate  that  fluid-mud 
der.  .i+y  undergoes  changes  of  up  to  0.05  gm/cm’ during  a wave  period,  approximately  4.5 
percent  of  overall  near-bottom  density;  to  produce  these  periodic  density  changes  requires 
the  suspension  and  deposition  of  a sediment  cloud  with  a concentration  of  50,000  mg/1. 
Inasmuch  as  present  theories  and  observations  cannot  explain  fully  the  rapid  settling  of  such  a 
sediment  cloud,  it  is  suggested  that  a collapsing  interface  is  formed  where  interparticle 
contact  is  maintained.  Additional  studies  must  be  undertaken  to  determine  the  extent  of 
vertical  oscillations  similar  to  those  found  in  Mississippi  Delta  muds. 

It  appears  that  mudbanks  undergo  substantial  change  in  time  and  space  from  the  longer 
term  suspension  and  redeposition  of  up  to  0.8  m of  fluid  mud  that  is  related  to  stage  of  the 
tide.  The  overall  high  suspended-sediment  concentrations  at  LW,  approximately  4000  mg/1  in 
coastal  surface  waters,  are  explained  by  this  suspension.  The  exchange  between  fluid  mud 
and  suspended  sediment  may  cause  a near-bottom  fluid  mud  density  decrease  of  from  1.16  to 
1.05  gm/cm^  in  a 3-hour  period  of  time. 
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The  presence  of  solitary  waves,  high  suspended-sediment  concentrations,  and  a rela- 
tively uniform  angle  of  wave  approach  throughout  the  year  may  lead  to  extraordinarily  high 
sediment  transport  rates.  If  waves  are  assumed  to  have  a net  drift  given  by  solitary  wave 
theory,  then  in  taking  angle  of  wave  approach  and  actual  measured  concentration  of  sediment 
in  suspension,  enough  sediment  can  be  transported  by  waves  alone  to  explain  the  propagation 
of  mudbanks  by  fluid-mud  transport.  Based  on  reasonable  wave  and  sediment  conditions,  3-70 
X 10*  m^  of  mud  can  be  transported  alongshore  each  year.  An  important  point  concerning 
high-volume  transport  rates  found  by  taking  this  approach  is  that  they  do  not  require  breaking 
waves  and  the  concept  of  radiation  stress  and  a nearshore  circulation  cell,  as  normally  applied 
to  sandy  coasts. 

Although  data  were  collected  on  one  mudbank  only,  on  the  basis  of  extensive  aerial 
reconnaissance  it  is  concluded  that  solitary  waves,  high  suspended-sediment  concentrations, 
and  similar  wave-mud  interactions  occur  along  other  parts  of  the  Surinam  coast  where  fluid 
mud  is  present  and  occur  perhaps  as  far  northwest  as  the  Orinoco  River  and  as  far  southeast 
as  the  Amazon  River.  A closer  examination  of  wave-mud  interactions  offers  possibilities  for 
substantially  increasing  and  quantifying  our  knowledge  of  sediment  transport  in  the  nearshore. 
Application  of  these  results  to  other  muddy  coasts  where  fluid  mud  is  present  and  suspended- 
sediment  concentrations  are  high  is  considered  appropriate. 
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APPENDIX  1 


Wave/Fluid-Mud  Pressure  Gage 


Fundamental  to  an  understanding  of  fluid-mud  dynamics  is  knowledge  of  the  details  of 
suspension  of  fluid  mud  by  waves.  Unsteadiness  of  the  sediment-suspending  fluid  precludes  the 
use  of  sampling  devices  that  are  designed  for  uniform,  steady  flows.  Therefore,  a continuously 
recording  prototype  instrument,  referred  to  here  as  a wave/fluid-mud  pressure-sensing  device, 
was  designed  and  constructed  by  the  technical  staff  at  Coastal  Studies  Institute.  This  instru- 
ment has  the  important  feature  of  being  capable  of  monitoring  remotely  both  mean  and 
fluctuating  components  of  near-bottom  or  subbottom  density  changes  and  at  the  same  time 
obtaining  a time  history  of  surface-water  waves  or  tidal  elevation  changes.  Although  pressure 
transducers  have  been  used  in  a variety  of  oceanographic  instruments  for  many  years,  such  a 
system  as  this  could  not  be  obtained  commercially. 

The  system  contains  three  basic  components:  two  Statham  model  PA-506  Amplibridge 
pressure  transducers;  an  electronic  readout  package;  and  a two-channel  Gould  Brush  strip- 
chart  recorder  (Fig.  I. la).  Each  pressure  transducer  has  a pressure  range  of  8.9632  x 10 '’to 
2.2753  X I O*"  dynes/cm''  (13-33  psia),  which  is  amplified  to  give  a 0-  to  5-volt  nominal  output 
over  this  range.  A 100-m  conductor  cable  connects  each  transducer  to  the  electronic  readout, 
where  pressure  fluctuations  app)ear  on  a panel  meter,  if  desired,  the  signal  can  be  recorded 
continuously  on  strip-chart  recorder.  Three  option  switches  are  available  on  the  readout:  a 
filter  selector  for  real-time  or  filtered  output;  a mode  selector  that  gives  the  option  of 
recording  waves  only,  waves  and  differential  density,  or  water  elevation  and  differential 
density;  and  a multiplier  selector  to  increase  resolution  when  taking  filtered  data  with  small 
density  changes.  The  readout  and  strip-chart  recorder  are  battery  powered.  Chart  speed 
range  is  from  1 mm/sec  to  50  mm/sec. 

The  concept  of  pressure  variation  in  a moving  fluid  is  the  basis  for  this  system. 
Pressure  under  a wave,  less  atmospheric  pressure,  is  given  by 


P cosh  [2Ti(z+h)/L]  H 

" PS— oYhTIihTL) — 2 


/2Trx 

V l'" 


The  first  term  in  this  equation  is  a dynamic  component  resulting  from  kinetic  energy  of  water 
motion;  the  second  term  is  a static  component  resulting  from  water  level  changes  associated 
with  the  wave  form.  Since  in  the  linearized  case  n = H/2  cos  (2itx/L  - 2vt/T),  equation  1.1  can 
be  written 

„ . _ cosh  [ 2ii  (z+h)/L  ] „ 
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This  can  now  be  expressed  as 


P = pg  (nK, 
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Figure  I.  la.  Wave/fluid-mud  pressure-sensing  device 
(recorder  not  shown). 
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Figure  I.lb.  Schematic  of  field  setup  for  wave/fluid-mud 
pressure-sensing  system  (see  Fig.  2). 


where  = [cosh  [2TT(z+h)/L]  ]/[cosh  (Zirh/L)}  is  the  pressure  response  factor.  In  shallow 
water  h/L  < 0.05,  K^-*- 1,  and  equation  1.3  reduces  to 

P (n- z)  =pgH^  . (1.4) 

This  is  identically  the  hydrostatic  pressure  equation  so  that,  when  dealing  with  solitary  waves, 
pressure  to  a first  approximation  is  that  of  the  static  component  only.  Values  of  at  the 

bottom  as  a function  of  water  depth  and  wave  length  are  tabulated  in  CERC  (1973)  and  show 
that  wave  attenuation  with  depth  in  shallow-water  situations  is  negligible.  Thus  the  difference 

in  hydrostatic  pressure  at  two  depths  below  the  surface,  H . and  Ho,  where  P,  =pgPa  and  Po 

Ad  a a d 

= pgHg,  can  be  written  as 


Pq  - = Pg(Hg  - H^)  . (1.5) 

If  Hg  - is  constant  with  time,  then  any  change  in  P^  - P^  is  directly  proportional  to  a 
density  change. 

Applying  the  above  hydrostatic  pressure  concept  in  the  field,  two  pressure  transducers, 
A and  B,  are  mounted  securely  to  a vertical  pipe  with  a known  distance  between  them.  Ho  - 

D 

(Fig.  I. lb).  An  adjustable  sleeve  arrangement  allows  the  depth  of  measurement  to  be  varied 

so  that,  if  desired,  both  transducers  can  be  raised  out  of  the  fluid  mud  to  record  wave  data 
only.  Hydrostatic  pressure  resulting  from  water  (or  fluid  mud)  head  is  nulled  out  of  each  pres- 
sure transducer  by  a potentiometer  located  on  the  instrument  panel.  At  the  beginning  of  each 
experiment,  differential  pressure,  P^  -P^j  zero;  only  a change  in  density  between  the 

pressure  transducers  will  cause  the  output  to  deviate  from  zero. 

Figure  2 shows  the  two  arrangements  of  pressure  transducers  relative  to  the  fluid-mud 
surface.  The  first  arrangement  is  used  to  monitor  surface  waves  on  one  channel  and 
simultaneous  density  changes  at  the  fluid-mud/water  interface  on  the  other  channel.  This  was 
the  normal  method  of  operation  at  site  1.  The  second  arrangement  is  used  to  record  long-term 
density  changes  in  fluid  mud  resulting  from  tidal  elevation  changes.  In  this  case  both 
transducers  are  placed  in  the  fluid  mud,  as  shown  at  site  2.  As  long  as  surface  waves  are 
absent  on  the  fluid-mud  surface,  the  upper  transducer  will  record  SWL  elevation  changes. 
When  used  in  either  fashion,  channel  1 is  recorded  in  real-time  (equation  1.4)  and  channel  2 is 
recorded  as  the  difference  between  filtered  A and  filtered  B (equation  1.5).  Best  results  were 
obtained  on  channel  2 by  using  a filter  with  20-sec  time  constant  and  a resolution  increase  of 
lOx. 


Prior  to  the  field  experiments  in  Surinam,  each  pressure  transducer  was  calibrated 
electronically  in  the  laboratory  so  that  an  input  voltage  equivalent  to  a force  of  10.05215  x lO"* 
dynes/cm^  (1.458  psi)  was  recorded  as  a 1-m  hydraulic  head.  This  is  the  force  exerted  by  a 
cubic  meter  of  clear  seawater  with  Te  = 15*0,  S = 33.70  7oo,  and  weight  density  = 1.025 
gm/cm^(64.0  Ib/ft^.  These  may  be  considered  "average"  seawater  conditions;  adjustments  in 
calibration  can  be  made  for  specific  localities. 

The  introduction  of  high  concentrations  of  sediment  into  seawater  will  increase  density. 
To  evaluate  the  effect  of  suspended  sediment  on  wave  height  as  determined  by  pressure- 
sensing devices  (equation  1.4),  Table  I.l  was  prepared  for  concentrations  of  a range  normally 
observed  in  Surinam.  Seawater  density  was  calculated  using  Te  = 30*0,  S = 34.75 ‘“/oo,  and 
specific  gravity  of  sediment  particles  = 2.65.  At  concentrations  less  than  50,000  mg/1,  the 
effect  on  density  is  small  enough  that  the  difference  between  apparent  and  actual  wave  height 
can  be  neglected.  Concentrations  greater  than  100,000  mg/1  may  significantly  affect  recorded 
wave  height;  however,  at  no  time  during  the  field  study  did  the  average  concentration  of  the 
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entire  water  column  change  by  100,000  mg/1  with  the  passage  of  waves. 


The  wave/fluid-mud  pressure  gage  was  field  calibrated  by  taking  samples  of  fluid  mud 
between  the  sensing  diaphrams  of  pressure  transducers  A and  B and  returning  them  to  the 
laboratory  for  bulk  density  determinations  (Appendix  III).  By  doing  this,  initial  bulk  density  can 
be  established  when  Pg  - = 0.  Any  long-term  deviations  from  the  zero  pen  position  are 

recorded  as  equivalent  water  height  changes  and  can  readily  be  converted  to  density  changes 
using  calibration  information  and  equation  1.5. 

Complications  arise,  however,  when  trying  to  determine  density  change  from  equivalent 
water  height  change  on  a short-term  wave-by-wave  basis.  The  RC  response  curve  for  a filter 
with  a 20-sec  time  constant  (Fig.  1. 2a)  shows  that  output  amplitude  for  waves  of  period  less 
than  five  time  constants  (100  sec)  is  less  than  100  percent  of  input  amplitude.  Thus  the  signal 
from  the  higher  frequency  waves  tends  to  be  attenuated.  Further,  output  amplitude  depends 
on  wave  shape,  so  that  in  natural  situations,  where  the  shape  of  each  wave  may  differ  from  the 
preceding  wave,  an  exact  differential  density  scale  cannot  be  determined.  For  waves  over 
mudbanks  in  Surinam,  the  scale  for  near-bottom  density  changes  resulting  from  passing  gravity 
waves  was  simplified  by  taking  an  idealized  cosine  wave  2 sec  in  duration  and  15  sec  in  crest- 
to-crest  period  (Fig.  1.2b), 


a.  RC  RESPONSE  CURVE 


3 


b.  INPUT  WAVES  FOR  DIFFERENTIAL  DENSITY  DETERMINATION 
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Figure  I.2a-b.  RC  response  curve  for  20-sec  time-constant  filter  (1. 2a);  input 
waves  used  for  output  amplitude/input  amplitude  ratio  determination  (1.2b). 


46 


Table  I.l 


Effect  of  Suspended  Sediment  on  Recorded  Wave  Height 
As  Determined  by  Pressure-Sensing  Instruments 


Sediment 
Concentration 
(10  mg/l) 

Density 
(gm/cm  ^ 

Error  in 

Recorded  Height 
(%) 

1 

1.02176 

0.3 

5 

1.02201 

0.3 

10 

1.02231 

0.3 

20 

1.02293 

0.2 

30 

1.02354 

0.1 

1.02416 

0.1 

50 

1.02477 

0.1 

100 

1.02784 

0.3 

200 

1.03415 

0.9 

300 

1.04038 

1.5 

400 

1.04038 

2.1 

500 

1.05283 

2.7 

1000 

1.08314 

5.3 

Precision  of  this  wave/fluid-mud  pressure  sensing  system  is  excellent.  According  to  the 
manufacturer,  total  error  in  pressure  transducers,  including  linearity  and  hysteresis,  is  less 
than  0.35  percent  over  the  13-33  psi  range.  The  electronic  readout  and  strip-chart  recorder 
each  have  a precision  of  1 percent.  If  maximum  error  for  the  system  is  taken  to  be  the  sum  of 
error  for  the  components,  then  results  are  reproducible  to  within  less  than  2.5  percent. 

Absolute  accuracy  of  the  system  is  unknown  but  is  believed  to  be  within  5 percent. 
Systematic  errors  in  differential  density  measurements  can  result  if  is  not  accurately 

determined.  Density  changes  resulting  from  the  passage  of  waves  are  believed  to  be  accurate 
to  within  0.008  gm/cm’,  while  long-term  density  changes  are  believed  to  be  accurate  to  within 
0.00^  gm/cm  \ 


APPENDIX  II 


Data  Processing  and  Time  Series  Analysis 


Digitizing 

Wave  and  differential  pressure  data  were  collected  as  analogue  signals  on  a strip-chart 
recorder.  Each  time  series,  of  20-min  record  length,  was  digitized  using  a Calmagraphic  III 
digitizer  (Calma  Corporation),  and  this  information  was  then  stored  on  magnetic  tape.  The 
digital  sampling  increment  of  At  = 0.2  sec  gave  6000  data  points  per  record  and  assured  a 
sufficiently  high  Nyquist  frequency,  f = l/2At,to  preclude  aliasing  effects  (Kinsman,  1965,  p. 
>*52).  Error  introduced  by  high-frequency  unsteady  hand  movements  was  assumed  negligible. 


Wave  Spectral  Analysis 

Spectral  analysis  is  an  analytical  tool  for  obtaining  essential  information  from  time 
series  records.  The  pioneering  work  by  Tukey  (1949,  1958,  1961)  formed  a base  from  which 
later  applications  were  taken  (Bendat  and  Piersol,  1966;  Kisel,  1969;  Edge  and  Liu,  1970). 

Given  a stationary  process  with  n (t)  as  the  randomly  fluctuating  signal,  then  the 
convolution  of  ri(t)  with  itself  at  a given  lag  t forms  the  autocorrelation  function,  defined  as 

, T 

r^(T)  / ri(t)r|(t+T)  dt  . (II.l) 

The  Fourier  transform  of  the  autocorrelation  function  is  the  sample  spectral  density  function 


R.j.fa)  = 2 / r^fr)  exp  (i2iTOT)  dx  . (11.2) 

o 

Variance,  defined  as  the  average  squared  deviation  of  observations  from  their  mean,  is  for  a 
process  of  zero  mean,  equal  to  the  average  Ti(t)^  Further, 

00 

Var  = J R^tr)  da.  (11.3) 

o 

Thus  total  variance  of  an  energy  spectrum  is  equal  to  the  total  area  under  all  spectral  peaks;  in 
small-amplitude  wave  theory,  this  energy  is  proportional  to  the  square  of  wave  amplitude. 
From  this  knowledge  it  is  possible  to  gain  important  information  about  the  distribution  of  wave 
energy  as  a function  of  frequency. 

Wave  spectral  analyses  were  performed  using  a Biomedical  Time  Series  Program 
(BMD02T),  modified  to  handle  ICT  data  points.  As  suggested  by  Kinsman  (1965,  p.  446),  the 
number  of  lags  selected,  m,  was  always  less  than  O.IN,  one-tenth  the  number  of  sample  points. 
Prior  to  analysis,  each  record  was  detrended  to  remove  the  effect  of  an  increase  or  decrease  in 
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water  depth  owing  to  changing  tide. 


Car*-  must  be  taken  when  interpreting  a wave  spectrum.  Although  an  underlying 
assumptic  is  stationarity,  in  reality  a 20-min  time  series  may  be  too  short  to  adequately 
characterize  a process  of  time  length  infinity.  Further,  difficulties  may  arise  when  trying  to 
resolve  sinusoidal  components  from  signals  that  are  strongly  nonsinusoidal.  With  periodic  sea 
surface  elevations  that  are  above  but  not  below  a SWL  "trough"  elevation,  as  in  pure  solitary 
waves,  harmonic  spectral  peaks  arise.  Although  the  variance  of  a wave  spectrum  is  a 
conservative  property,  the  derived  height  information  is  not,  and  harmonic  energy  peaks  do  not 
represent  true  sea  surface  elevations. 


APPENDIX  III 


Sample  Analysis 


Water  Samples 

Water  samples  were  collected  in  a 1-liter  modified  Van  Dorn  bottle,  stored  in  airtight 
500-ml  plastic  containers,  then  returned  to  the  laboratory  for  analysis.  Suspended -sediment 
concentration  and  salinity  were  determined  for  each  sample. 

Samples  of  30-300  ml  were  filtered  through  preweighed  0.45micron,  47-mm  diameter 
Millipore  membrane  filters  using  a pressure  filtration  system.  Preliminary  laboratory 
investigations  indicated  that  a stacked  filter  system  to  account  for  leaching  was  not  necessary 
owing  to  the  small  volumes  of  water  filtered  (Eaton  et  al.,  1969).  After  filtering,  each  filter 
pad  was  washed  three  times  with  10  ml  of  distilled  water  to  remove  salts,  air  dried  at  25®C  for 
48  hours,  desiccated  for  1 hour,  and  weighed  to  the  nearest  0.01  mg.  Suspended  sediment 
concentration  was  then  calculated  and  expressed  as  mass  of  dried  sediment  per  volume  of 
solution.  The  precision  of  this  method  depends  on  concentration  of  the  sample;  however,  the 
coefficient  of  variation  (standard  deviation/mean)  was  less  than  3 percent. 

Salinity  was  determined  from  the  filtrate  of  each  sample  using  a Mark  IV  conductive 
salinometer  (Lab  Line  Instruments,  Inc.).  With  this  instrument,  accuracy  is  a function  of 
conductivity.  According  to  the  manufacturer,  accuracy  is  better  than  -0.5  "/»  at  35  “/oo  and 
-O.l  “oo/  at  10  7oo  . 


Fluid-Mud  Samples 

"Suspensions"  with  concentrations  greater  than  10,000  mg/1  are  generally  referred  to  as 
fluid  mud  (Krone,  1962)  and  cannot  be  analyzed  by  standard  filtration  techniques.  To 
determine  concentration  of  fluid-mud  samples,  50  ml  of  solution  were  carefully  measured  from 
a volumetric  flask,  rinsed  four  times  witn  distilled  water  to  remove  salts,  the  excess  being 
decanted  each  time,  then  oven  dried  at  80“  C.  After  cooling,  the  sample  was  weighed  at  room 
humidity  to  the  nearest  0.01  mg  and  the  concentration  expressed  as  mass  of  solids  per  volume 
of  solution. 

Bulk  density,  defined  as  total  mass  of  sample  per  volume  of  solution,  can  be  determined 
in  the  laboratory  by  weighing  a known  volume  of  fluid  mud.  However,  precision  is  f>oor  since 
small  volume  errors  can  introduce  large  weight  errors.  Overall  precision  is  improved  by  taking 
the  specific  volume  of  sediment  (0.377  cmVgm  for  particle  density  of  2.65  gm/cm’)  and 
determining  the  amount  of  seawater  that  will  be  displaced  for  a given  concentration,  as  deter- 
mined above.  For  example,  sediment  particles  in  fluid  mud  with  a concentration  of  100,000 
mg/1  in  seawater  of  density  1.025  gm/cm*  will  displace  38.6783  gm  of  water,  assuming  particle 
density  is  2.65  gm/cm’.  Bulk  density  is  then  calculated  (per  liter  of  solution)  as 
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(initial  mass  of  seawater)  - (mass  of  water  displaced) 

1025  gm  - 38.68  gm 

+ (mass  of  sediment) 

100  gm  = 1086  gm 

This  corresponds  to  a bulk  density  of  1.086  gm/cm\  Following  a similar  procedure,  water 
content,  given  as 


[(mass  of  sample  - mass  of  solids)/(mass  of  solids  + salts)  ] x 100 
[ (1086  gm  - 100  gm)  ^(100  gm  + 25  gm)]  x 100 
= 709% 


(III. 2) 


can  be  determined  without  weighing  a measured  volume.  Mass  of  salts  must  be  known  from  a 
salinity  determination.  Since  the  void  spaces  in  a marine  sediment  are  assumed  to  be  water- 
saturated,  voids  ratio  is  given  by 


Volume  of  voids 
Volume  of  solids 


Specific  gravity  x water  content 
100 


(III.3) 


Richards  et  al.  (197^)  have  found,  from  analysis  of  more  than  1500  marine  sediment  samples, 
that  the  linear  relationship  between  voids  ratio  and  water  content  gives  an  average  specific 
gravity  of  2.72.  ^ 

Viscosity  was  determined  graphically  from  curves  prepared  by  NEDECO  (1965)  that 
were  based  on  more  than  UOOO  viscosity  measurements  made  on  fluid  mud  and  soil  samples 
from  Surinam  and  British  Guiana  using  a Brookfield  viscometer.  Viscosity  is  a function  of 
sediment  concentration,  salinity,  and  shearing  rate.  Values  reported  in  the  present  study  are 
for  seawater  of  salinity  35  7oo  and  shearing  rate  of  10  rpm  for  a duration  of  6 min. 

Particle  size  distributions  were  determined  from  pipette  analyses  following  standard 
procedures  given  by  Folk  (1968).  Approximately  15  gm  of  sediment  (wet  weight)  were 
disaggregated  in  a blender  filled  halfway  with  10  percent  calgon  dispersant  solution,  then 
poured  into  a graduated  cylinder  and  brought  to  exactly  the  1000- 
ml  mark  with  dispersant  solution.  Each  cylinder  was  left  for  24  hours  to  check  for 
flocculation.  Withdrawals  were  made  at  appropriate  times  and  depths  (assuming  Stokes 
settling  velocity)  and  released  into  evaporating  dishes  for  drying.  Pipette  withdrawals  were 
terminated  at  1 l^i  (0.5  micron),  since  a constant  temperature  could  not  be  assured  past  this 
point.  Each  sample  was  desiccated  and  weighed,  a calgon  correction  was  subtracted,  and  the 
weight  percentage  was  determined. 
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Time  senes  measurements  of  shallow-water  waves  and  fluid-mud  density  variations,  taken 
simultaneously  with  tide  elevation  and  suspended  sediment  data,  indicate  that  wave/fiuid-mud 
interactions  in  'lie  nearshore  may  be  largely  responsible  for  the  present-day  accumulation  of  fine- 
grained sediment  on  the  open,  unprotected  coast  of  northeastern  South  America.  Results  ol  field 
experiments  conducted  along  the  central  Surinam  coast  show  that  accumulations  of  fluid  mud 
affect  incoming  swell  by  changing  their  form  from  sinusoidal  to  solitary-like  and  by  preventing 
wave  breaking  except  for  occasional  spilhng.  As  long'Penod  swell  (T  « 12-16  sec)  propagates  over 
shallow  banks  of  fluid  mud,  the  wave  height  to  water  depth  ratio  remains  nearly  constant  at  0.23; 
the  steady  decrease  in  wave  height  with  shoaling  water  depth  indicates  that  substantial  amounts  of 
wave  energy  are  lost  to  a fluid-mud  bottom  even  though  breaking  does  not  occur.  Evaluation  of 
several  mechanisms  for  the  dissipation  of  wave  energy  reveals  that  attenuation  of  high-frequency 
waves  (T  < 3 sec)  may  be  due  to  internal  friction  resulting  from  the  extraordinarily  high  kinematic 
viscosity  of  muddy  coastal  waters,  whereas  longer  period  waves  may  lose  energy  to  viscous  shear 
within  a bottom  boundary  layer.  (U) 

A specially  designed  wave/lluid-mud  pressure-sensing  system  indicates  that  fluid-mud 
density  fluctuates  by  0.023  to  0.0)  gm/cm*  on  a wave-by-wave  basis.  As  each  wave  crest  passes 
the  sensing  instrument,  a cloud  of  sediment  with  particle  concentration  of  approximately  50,000 
mg/1  is  suspended  instantly;  rapid  settling  that  follows  is  explained  by  the  formation  of  a collapsing 
interface.  Wave-by-wave  sediment  suspension  is  superimposed  on  a longer  term  suspension  and 
deposition  process  that  is  related  to  stage  of  the  tide.  Surface  suspensate  concentrations  may 
exceed  3500  mg/I  during  this  exchange  process.  (U) 

Utilizing  the  concept  of  wave-associated  currents,  it  has  been  shown,  in  taking  angle  of 
wave  approach,  average  observed  concentration  of  sediment  in  suspension,  and  net  drift  as  given  by 
solitary  wave  theory,  that  volume  transports  of  3 to  70  x 10*  mVyr  can  be  explained  by  waves 
alone.  If  a hypothesis  of  mud  transport  by  solitary  waves  is  accepted,  then  the  continuous 
shoreward  transport  of  suspended  fluid  mud,  combined  with  the  high  rate  of  wave  energy 
dissipation,  explains  in  part  the  ability  of  muddy  coasts  to  protect  their  shorelines.  (U) 
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